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In this paper we consider the simplified form that a recently introduced general operator descrip- 
tion of the Hubbard model on the square lattice with N% S> 1 sites, effective transfer integral t, and 
onsite repulsion U has in a suitable one- and two-electron subspace. Such an operator description is 
that consistent with the model exact global symmetry recently extended to SO(3) x SO(3) x U(l). 
Our study profits from the rotated-electron occupancy configurations in the original lattice that 
generate the energy eigenstates being described by occupancy configurations of spin- 1/2 spinons in 
a spin effective lattice, r;-spin-l/2 77-spinons in a ry-spin effective lattice, and c fermions in a c effec- 
tive lattice. The latter three types of occupancy configurations refer to state representations of the 
spin SU(2) symmetry, jy-spin SU(2) symmetry, and hidden U(l) symmetry recently found in Ref. 4 , 
respectively. In the limit of very large number of lattice sites N% S> 1 that such a description refers 
to the emergence of the above three independent effective lattices simplifies the study of the effects 
of hole doping on the spin subsystem. In the one- and two-electron subspace the model refers to a 
square-lattice quantum liquid that plays the same role for the Hubbard model on a square lattice as 
the Fermi liquid for isotropic three-dimensional correlated perturbative models. There is a large con- 
sensus that in the thermodynamic limit N% — ¥ oo long-range antiferromagnetic order occurs in the 
spin-density m = ground state of the half-filled Hubbard model on the square lattice. Here we find 
that the corresponding spontaneous symmetry breaking lowers the symmetry of such a state from 
SO(3) x £0(3) x 17(1) for 7V a 2 > 1 large but finite to [17(2) x U(l)]/Z$ = [SO (3) x (7(1) x U(l)]/Z 2 
for —¥ oo. Moreover, we argue that the spin effective lattice being identical to the original lat- 
tice is a necessary condition for the occurrence of ground-state long-range antiferromagnetic order 
in the limit — > oo. Consistently, strong evidence is provided that for very small hole concen- 
tration < x <C 1 the ground state has a short-range incommensurate-spiral spin order. (The 
related investigations of Ref.— provide evidence that a spin short-range order exists for < x < x», 
whereas for x > x, the ground state is a spin disordered state. Here i, > 0.23 for approximately 
U /At > 1.3.) Our results are of interest both for condensed-matter systems and ultra-cold fermionic 
atoms on an optical square lattice. Elsewhere evidence is provided that upon addition of a weak 
three-dimensional uniaxial anisotropy perturbation to the square-lattice quantum liquid, its short- 
range spin order coexists for iV„ — ¥ oo, low temperatures, and a well-defined range of finite hole 
concentrations with a long-range superconducting order. 

PACS numbers: 71.10.Fd, 71.10.-fw, 71.27. +a, 71.10.Hf, 71.30.+h 

I. INTRODUCTION 

The Hubbard model on a square lattice is the simplest toy model for describing the effects of electronics correla- 
tions in the high-T c superconductors^ and their Mott-Hubbard insulators parent compounds-. In the case of such 
superconductors, addition of a weak three-dimensional (3D) uniaxial anisotropy perturbation associated with a small 
effective transverse transfer integral t± is required to capture some of the basic properties of their physics. 

The studies of this paper focus on the model on a square lattice. It has no exact solution, so that many open 
questions about its properties remain unsolved. A recent exact result, which applies to the model on any bipartite 
lattice!, is that for on-site repulsion U > its global symmetry is SO (3) x SO(3) x U(l) = [SO {A) x U(l)]/Z 2 . That 
is an extension of the model well-known SO (4) symmetry^, which becomes explicit provided that one describes the 
problem in terms of the rotated electrons obtained from the electrons by any of the unitary transformations of the 
type considered in Refs^i Another exact result useful for the studies of this paper is that the ground state of the 
repulsive half-filled Hubbard model on a bipartite lattice is for finite number of lattice sites and spin density m = 
a spin singled. 

The model can be experimentally realized with unprecedented precision in systems of ultra-cold fermionic atoms on 
an optical lattice with any geometry^. In the experimental context of cold-fermion optical lattices, Ref^ discusses 
the possibilities to approach the pseudogap or ordered phases by manipulating the scattering length or the strength 
of the laser-induced lattice potential. In turn, the large square lattices considered in RefJ^ for the half-filled model 
provide improved resolution of the Greens function in momentum space, allowing a more quantitative comparison with 
time-of-flight optical lattice experiments. Such large scale determinant quantum Monte Carlo (DQMC) calculations 
provide as well useful information on the effective bandwidth, momentum distribution, and magnetic correlations of 
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the half-lling model. 

A. The model 

The studies of this paper profit from the simplified form that the general description of the Hubbard model on 
the square lattice with N% = [N a ] 2 ^> 1 sites, effective nearest- neighbor transfer integral t, and onsite repulsion U 
introduced in Rcf.i 2 - has in a suitable one- and two-electron subspace considered below. The hole concentration and 
spin density read x = [N% — N]/N% and m = [N^ — N$/N%, respectively. The model on such a lattice with spacing 
a, N a 3> 1 even, and lattice edge length L = N a a is given by, 

N l 

H = tf + U[Nl-Q]/2 ] f = -Y, + Q = EE ^(l-nfj.-a). (1) 

Here T is the kinetic-energy operator in units of t and Q is the operator that counts the number of electron singly 
occupied sites so that the operator D — [N — Q]/2 counts the number of electron doubly occupied sites, n? j a = 

c f-,a c ?i,° where ~ a =t (and -a =i) for a =1 (and a =f), AT = N <?i and ^ = J2j=i n ? 3 ,<?- 
The kinetic-energy operator T given in Eq. ([1} can be expressed in terms of the operators, 

T = - E Yl{ rif 'j-° c ? j ,<T C f' j >," n f' j >-* + i 1 - npj-^ct^cp.,^ (1 - ?V 3 ,,- CT )] , 

as T = To + +T-x- As discussed in Refs.— i 12 i 13 , these three kinetic operators play an important role in the physics. 
The operator To does not change electron double occupancy whereas the operators T + \ and T_i do it by +1 and —1, 
respectively. 

B. The general rotated-electron operator description 

The general operator representation introduced in RefJ^ applies for N% 3> 1 and refers to suitable quantum objects 
whose occupancy configurations generate the state representations of the group 5*0(3) x 5*0(3) x U(l). Addition of 
chemical-potential and magnetic-field operator terms to the Hamiltonian lowers the model symmetry. Such terms 
commute with it so that its 4 N °- momentum and energy eigenstates correspond to representations of that group for all 
values of x and m. It is shown in such a reference that all the physics of the model in the whole Hilbert space can be 
obtained from that of the model in the subspace spanned by the lowest-weight states (LWSs) of both the 77-spin SU (2) 
and spin SU (2) algebras. In this paper we profit from the simpler properties that the general operator description 
introduced in Rcf. 12 has in a suitable one- and two- electron subspace. That subspace is spanned by a well-defined set 
of energy eigenstates, which for U/At > are generated by simple momentum occupancy configurations of spin-less 
charge objects and spin-singlet objects, respectively. 

The general description introduced in Ref refers for U /At > to a particular choice of the complete set of 4 N " 
energy, momentum, r;-spin, 7;-spin projection, spin, and spin-projection eigenstates {\^u/At)}- I n the limit U/4t — > 00 
such states correspond to one of the many choices of sets {l^oo)} of 4 W » U / At-inSmte energy eigenstates. For the 
choice corresponding to the description of Refji 2 - there exists exactly one unitary operator V — V(U /At) such that 
\^u/4t) = ^^|*oo)- The point is that for most choices of the set of energy eigenstates {l^oo)} the corresponding states 
Y&U/At) = ^^1*00) are not energy eigenstates for U/At > 0. For the description introduced in that reference such 
states are energy eigenstates. The corresponding unitary operator V = V(U/4t) and the rotated-electron operators, 

4j ,a = ^ 4j ,(r V ; c F . , CT = t> f c F . , CT V ; n P] , ff = ct a c P] >CT , (3) 

play a central role in that description. The states \^u/At) — ^l^oo) (one for each value of U/At > 0) that are 
generated from the same initial state \^oo) belong to the same V tower. 
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Any operator can be written in terms of rotated-electron creation and annihilation operators as, 

6 = vov* = o + [d, s] + ^[[d, s], §} + ... = 6 + [d, §} + ^[[d, §}, §} + ..., 



t 



0{t 2 /U z ); S = -- T + i-T_! +0(t*/U J ). (4) 



Here the operator O = V' OV has the same expression in terms of rotated-electron creation and annihilation operators 
as O in terms of electron creation and annihilation operators, respectively. The operator 5 appearing in Eq. Q is 
related to the unitary operator as V' = e s and V = e~ s . Since for finite U/4t values the Hamiltonian H of Eq. 
(fTJ) does not commute with the unitary operator V = e~ s , when expressed in terms of rotated-electron creation and 
annihilation operators it has an infinite number of terms. According to Eq. Q it then reads, 

H = VHV* =H+[H,S] + - [[H, S], S] + ... . (5) 

The commutator [H, 5] does not vanish except for U/At — > oo so that H ^ H for finite values of U/At. 

For U I At very large the Hamiltonian of Eq. ([5]) corresponds in terms of rotated-electron creation and annihilation 
operators to a simple rotated-electron t — J model. In turn, the higher-order t/U terms become increasingly important 
upon decreasing U /At. They generate effective rotated-electron hopping between second, third, and more distant 
neighboring sites. Indeed, the products of the kinetic operators To, T + \, and T_i contained in the higher-order terms 
of 5 = —(t/U) [T+i — T_i] + 0(t 2 /U 2 ) also appear in the Hamiltonian expression ([5]) in terms of rotated-electron 
creation and annihilation operators of Eq. ([3]). In spite of the operators Tq, T +1 , and T_i generating only rotated- 
electron hopping between nearest-neighboring sites, their products generate effective hopping between for instance 
second and third neighboring sites. For instance, the Hamiltonian terms generated up to fourth order in t/U are within 
a unitary transformation the equivalent to the t—J model with ring exchange and various correlated hoppings 13 . 
The real-space distance in units of the lattice spacing a associated with the effective hopping between second and 
third neighboring sites is for the model on the square lattice y/2a and 2 a, usually associated with transfer integrals 
t' and t", respectively. For instance, for the half filled Hubbard model on the square lattice in the subspace with 
both vanishing rotated-electron double occupancy and vanishing rotated-hole double occupancy the interactions can 
be expressed completely by spin operators. The corresponding Hamiltonian expression is given below in Section III-B 
up to fifth order in t/U. Some of the interactions in spin space in its t /U 3 terms refer indeed to pair of spins at 
second and third neighboring sites. 

It follows that when expressed in terms of the rotated-electron operators emerging from the specific unitary trans- 
formation considered above the simple Hubbard model (JT|) as given in Eq. (|5|) contains Hamiltonian terms associated 
with higher-order contributions. Those can be effectively described by transfer integrals t', t", and of higher order. 
For hole concentration equal to or larger than zero both the model ground state and the excited states that span the 
one- and two-electron subspace considered below have vanishing rotated-electron double occupancy. Fortunately, for 
intermediate and large values of U/At obeying approximately the inequality U/At > uq ~ 1.3 and thus t/U < 0.2, 
besides the original nearest-neighboring hopping processes only those involving second and third neighboring sites are 
relevant for the square-lattice quantum liquid. That quantum liquid refers to the Hamiltonian of Eqs. ([T]) and ([5]) 
in the one- and two-electron subspace. The value U/At — u w 1.302 is that at which an important energy scale A 
whose U / At dependence is studied below in Section III-C reaches its maximum magnitude. 

Hence for approximately U/At > uq, out of the infinite terms on the right- hand-side of Eq. (|5|) only the first few 
Hamiltonian terms play an active role in the physics of the Hubbard model on the square lattice in the one- and 
two-electron subspace. Therefore, for intermediate and large values of U/At such a square-lattice quantum liquid 
can be mapped onto an effective t — J model on a square lattice with t, t 1 — t'(U/At), and t" = t"(U/At) transfer 
integrals. The role of the processes associated with t 1 — t'(U/At) and t" = t"(U/At) becomes increasingly important 
upon decreasing the U/At value. 



C. Three basic objects emerging from the rotated-electron description 

In addition to the spin S s and ?7-spin associated with the spin SU(2) and 77-spin SU(2) symmetries, respectively, 
of the model global 50(3) x 50(3) x U(l) = [SU(2) x 577(2) x U(l)]/Z 2 symmetry, the eigenvalue 5 C of the generator 
of the hidden U(l) symmetry found in Refi^ plays an important role in the general rotated-electron description of 
RefJ^. It is such that 25 c is the number of rotated-electron singly occupied sites, which is a good quantum number 
for U/At > 0. We denote the 77-spin and spin projection by 5^ 3 = — [N 2 — N]/2 and 5f 3 = — [JVf — ATJ/2, respectively. 
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That the global symmetry of the Hubbard model on a bipartite lattice is [SU(2) x 5/7(2) x t/(l)]/Z| rather than 
SU(2) x SU(2) x U(l) is due to only state representations such that 2S V , 2S S , and 2S C are simultaneously even or 
odd integer numbers being allowed^. The group SU(2) x SU(2) x /7(1) has the same seven generators but four times 
more representations than the group 50(3) x 50(3) x U(l). 

Studies of the Hubbard model on the square lattice that rely on its transformation properties under symmetry 
operations 14 should be extended to account for the new hidden U(l) symmetry recently found in Ref.—. Indeed that 
the global symmetry of the Hubbard model on a square lattice and all other bipartite lattices is larger than SO (4) 
and given by SO(3) x SO (3) x /7(1) is expected to have important physical consequences. In the case of the one- 
dimensional (ID) bipartite lattice the model has an exact solutionis—. The main reason why its solution by the 
algebraic Bethe-ansatz inverse scattering method^ was achieved only thirty years after that of the coordinate Bethe 
ansat a 15 i 16 is that it was expected that the charge and spin monodromy matrices had the same traditional ABCD 
form, found previously for the related ID isotropic spin 1/2 Heinsenberg models. Such an expectation was that 
consistent with the occurrence of a spin SU(2) symmetry and a charge (and 77-spin) SU(2) symmetry known long 
ago£, associated with a global SO (A) — [SU(2) x SU(2)]/Z 2 symmetry. If that was the whole global symmetry of 
the ID Hubbard model, the charge and spin sectors would be associated with the 77-spin 5/7(2) symmetry and spin 
SU(2) symmetry, respectively. A global 50(4) = [5/7(2) x SU(2)]/Z 2 symmetry would then imply that the charge 
and spin monodromy matrices had indeed the same Faddeev-Zamolodchikov ABCD form. Fortunately, Martins and 
Ramos used an appropriate representation of the charge and spin monodromy matrices, which allows for possible 
hidden symmetries. For the particular case of the bipartite ID lattice the results of Refi^ reveal that the hidden 
symmetry beyond 50(4) is the charge global U{\) symmetry found in that reference. For U/At > the model charge 
and spin degrees of freedom are then associated with U(2) — SU(2) x /7(1) and 5/7(2) symmetries, rather than with 
two 5/7(2) symmetries, respectively. The occurrence of such charge U(2) — 5/7(2) x U(l) symmetry and spin 5/7(2) 
symmetry is behind the different ABCDF and ABCD forms of the charge and spin monodromy matrices of RefJ^, 
respectively. 

Addition of chemical-potential and magnetic-field operator terms to the Hamiltonian (fT]) lowers its symmetry. As 
mentioned above, such terms commute with it so that the global symmetry being [5/7(2) x 5/7(2) x U(l)]/Z 2 = 
50(3) x 50(3) x U(l) implies that a set of 4^ independent rotated-electron occupancy configurations generate the 
corresponding state representations of that global symmetry for all values of the electronic density n — (1 — x) and 
spin density m. Consistently, in Ref. 4 it is confirmed that the total number of such independent representations 
equals indeed the Hilbert-spacc dimension 4 N ° . The rotated electron occupancy configurations that generate such 
state representations are simpler to describe in terms of those of suitable related quantum objects. The eigenvalue 5 C 
of the new hidden U(l) symmetry controls the numbers of such objects, which are as well good quantum numbers. 
Indeed the investigations of Ref.— reveal the emergence within the rotated-electron description of three basic quantum 
objects: M s = 2S C spin-1/2 spinons, M n = [N% — 2S C ] 77-spin- 1/2 77-spinons, and N c — 2S C spin-less and 77-spin-less 
charge c fermions. The latter live on a lattice with N% = [N c + AT™] sites identical to the original lattice. Here 
N£ = [N% — 2S C ] gives the number of c fermion holes. 

The relation of such objects to the rotated electrons is as follows. The M s — 2S C spin-1/2 spinons describe the 
spin degrees of freedom of the 25 c rotated electrons that singly occupied sites. The M v = [N% — 2S C ] 77-spin-l/2 
77-spinons describe the 77-spin degrees of freedom of the [N% — 2S C ] sites doubly occupied and unoccupied by rotated 
electrons. Specifically, the 77-spinons of 77-spin projection —1/2 and +1/2 refer to the sites doubly occupied and 
unoccupied, respectively, by rotated electrons. The counting of the number of spinon and 77-spinon independent 
occupancy configurations is an exercise fulfilled in Ref A However, the internal structure of the specific spinon and 
77-spinon occupancy configurations of the 4 N ° U /At > energy eigenstates associated with the RefJ^ description is a 
complex problem, which partially simplifies in the one- and two-electron subspace considered below. In turn, in each 
rotated-electron configuration of a U /At > energy eigenstate, the iV c = 25 c c fermions exactly occupy the same 
sites as the 25 c rotated electrons that singly occupy sites. Furthermore, the N% = [N% — 2S C ] sites unoccupied by 
c fermions are those doubly occupied and unoccupied by the rotated electrons. The electronic charges of the 25 c 
rotated electrons that singly occupy sites are carried by the N c = 25 c c fermions whereas their spins are carried by 
the M s = 2S C spin-1/2 spinons. 

For the general description of RefJ^ the rotated-electron occupancy configurations referring to (i) the spin degrees 
of freedom of the singly occupied and (ii) the 77-spin degrees of freedom of the unoccupied and doubly-occupied sites 
are independent. They correspond to the state representations of the spin 5/7(2) symmetry M s = 2S C spin-1/2 
spinons and 77-spin 5/7(2) symmetry M v = 25' 1 77-spin-l/2 77-spinons, respectively. In turn, the U(l) symmetry state 
representations refer to the relative occupancy configurations of the 25 c rotated-electron singly-occupied sites and 25^ 
rotated-electron unoccupied and doubly-occupied sites. Those correspond to the c fermion occupancy configurations. 

Within the LWS representation of the general operator description introduced in RefJ^, the c fermion creation 
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operator has the following expression in terms of the rotated-electron operators of Eq. , 

/i,c = 4t( 1 -^^ + ei ^^.t^4; 4 

Here we have introduced the corresponding c fermion momentum-dependent operators as well and e" 1 '"^ is ±1 de- 
pending on which sub-lattice site fj is on. The use of Eq. ([3]) allows the expression of the c fermion operators ((6J in 
terms of electron operators. This involves the electron - rotated-electron unitary operator V. The c momentum band 
is studied in RefJ^ and has the same shape and momentum area as the electronic first-Brillouin zone. 

In turn, the three spinon local operators st and three 77-spinon local operators p 1 -. such that I — ±, z are given by, 

4 3 = nf^c q fj ; p l r s = (1 - n? jtC ) q^ , l = ±,x 3 . (7) 

Here q% 3 and = q% x ± i q% 2 are the rotated quasi-spin operators ql = si + nl . (We denote the Cartesian 
coordinates by x^,x 2 ,x 3 .) In addition, 

"* i ,c = f? jtD f? i ,c, (8) 

is the c fermion local density operator. In terms of rotated-electron creation and annihilation operators the rotated 
quasi-spin operators read, 

4= J, S+pJ = (4t-e i ^^,t)^; ^ = (4) f ; ^ = ^+^ = ^"^,1- (9) 

D. Two qualitatively different types of spinon (and 77-spinon) configurations and our notation 

An important result of RefJ^ is that a well-defined number of spin-1/2 spinons and n-spin-1/2 77-spinons remain 
invariant under the electron - rotated-electron unitary transformation. We call those independent spinons and in- 
dependent ?y-spinons, respectively. In the one- and two-electron subspace considered in this paper there are no 
77-spin-projection —1/2 77-spinons. The independent spinons and the independent 77-spin-projection +1/2 77-spinons 
have vanishing energy and momentum. There are no independent spinons in m = ground states. 

The spinons that are not invariant under the electron - rotated-electron unitary transformation are confined within 
spin-neutral 27/-spinon composite sv fermions. Here v = 1,2,... is the number of spinon pairs confined within a 
spin-neutral 2^-spinon composite sv fermion. The description of Ref used in our studies refers to a number of sites 
N~ 3> 1 very large but finite. Within such a description all m = ground-state spinons are confined within spin- 
neutral two-spinon si fermions. The condition that in an excited state spinons are deconfined is that they are invariant 
under the electron - rotated-electron unitary transformation. For instance, a spin-triplet excited state involves two 
independent spinons that correspond to an isolated mode below a continuum of two-sl-fermion-hole excitations. 

For a given state, the values of the numbers L v .±1/2 of independent ±1/2 77-spinons and L s ,±1/2 of independent 
±1/2 spinons are fully determined by those of the 77-spin S v and 77-spin projection S^ 3 = —x N% and spin S s and spin 
projection S* 3 — —mN%, respectively, as follows, 

L a = [L a ,-i/ 2 + L at+1 / 2 ] = 2S a ; L a ,±i/2 = [S a TS?]; a = ?7,s. (10) 

The invariance of such independent 77-spinons (and spinons) stems from the off diagonal generators of the 77-spin (and 
spin) algebra, which flip their 77-spin (and spin), commuting with the unitary operator V. Hence it follows from Eq. 
(HI that such generators have for [7/4/j > the same expressions in terms of electron and rotated-electron operators. 

Alike the spin-1/2 spinons, the 7y-spin-l/2 77-spinons that are not invariant under the electron - rotated-electron 
unitary transformation are confined within 77-spin-neutral 27/-f7-spinon composite 77^ fermions. Within the notation of 

Ref. 12 also used in this paper, v = 1, C n and v = 1 C s are the numbers of //-spinon and spinon pairs confined 

within a composite r\v fermion and a sv fermion, respectively. Their maximum values C v and C s , respectively, are 
given in Eq. (|A5|) of Appendix A. In principle there are some occupancy configurations of such objects that generate 
the exact momentum and energy eigenstates for U/At > 0, yet the detailed internal structure of such configurations 
is a complex problem^. 

Accordingly to the studies of RefJ^ that refer to N% ^ 1 very large but finite the m — and x > ground states 
are spin singlets with N s \ = M s /2 = S c = N/2 two-spinon si fermions, N c = 2S C = N = (1 — x)N% c fermions, 



N l 

= V e +,;9 > r > ft . (6) 

V a y = i 
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L v ,+i/2 = [N% — 2S C ] = x N% 77-spinon-projection +1/2 independent 77-spinons, no independent spinons, no sv fermions 
with v > 1 spinon pairs, and no r\v fermions. It is found below that for the Hubbard model on the square lattice 
in the one- and two-electron subspace considered in the studies of this paper only the c fermions and the two-spinon 
si fermions play an active role. In the two quantum liquids referring to (i) x = and m = and (ii) x > and 
rn = 0, respectively, considered in our studies all the ground-state M s = 2S C spinons are for N% ^> 1 very large but 
finite confined within N s i = M s /2 spin-neutral two-spinon si fermions. Spinon confinement and the corresponding 
deconfined degrees of freedom in second-order phase transitions is a problem of physical interest^ - — . 

In summary, the degrees of freedom of the rotated-electron occupancy configurations of each of the sets of N% = 

[N% — 2S C ] and N% = 2S C sites of the original lattice that generate the energy eigenstates separate into two types 
of configurations. A first type of occupancy configurations are those of the c fermions associated with the operators 
ft c of Eq. ([U where j = 1, ...,iVf. Such occupancy configurations correspond to the state representations of 

the global U(l) symmetry found in Ref4. The c fermions live on a c effective lattice, which is identical to the 
original lattice. The N c = N% = 2S C c fermions occupy the sites singly occupied by the rotated electrons. The 
— N% — [N% — 2S c ] rotated-electron doubly-occupied and unoccupied sites are those unoccupied by the c fermions. 
The c fermion occupancy configurations describe the relative positions in the original lattice of the N% = [N% — 25 c ] 
sites of the 77-spin effective lattice and N% = 2S C sites of the spin effective lattice. 

Consistently, the remaining degrees of freedom of rotated-electron occupancies of the sets of N% = [JV| — 2S C ] and 
N% b = 25 c original-lattice sites correspond to a second type of occupancy configurations. Those are associated with 
the 77-spin SU(2) symmetry and spin SU(2) symmetry representations, respectively. The occupancy configurations 
of the set of N% = [N% — 2S C ] sites of the 77-spin effective lattice and set of N% s — 2S C sites of the spin effective 
lattice are independent. The former configurations refer to the operators p^. of Eq. 0, which act only onto the 
N% = [N% — 2S C ] sites of the 77-spin effective lattice. The latter configurations correspond to the operators si. given 

in the same equation, which act onto the N% = 2S C sites of the spin effective lattice, respectively. This is assured by 
the operators (1 — np. )C ) and n? jtC in their expressions provided in that equation, which play the role of projectors 
onto the 77-spin and spin effective lattice, respectively. 

The spin-singlet two-spinon composite local si fermions of the description of RefJ^ live on their own si effective 
lattice. One occupied site of such a si effective lattice involves two sites of the spin effective lattice. Consistently, 
for spin density m — and spin S s = states of the one- and two-electron subspace considered below the number 
of sites of the si effective lattice reads N% = N% /2 = S c . As justified in Refs i 12 ' 31 , for N% — > 00 both the spin 
and si effective lattices are square lattices with spacing a s and a s \, respectively, related as a s iV2a s . Some further 
basic information on the general description introduced in Ref.— needed for the studies of this paper is provided in 
Appendix A. 

E. The goals and organization of this paper 

Our general aim is to show that in the one- and two-electron subspace considered below the physics of the Hubbard 
model on the square lattice simplifies and refers to the square-lattice quantum liquid further studied in Ref^. The 
set of energy eigenstates that span such a subspace are generated by momentum occupancy configurations of the c 
and si fermions. Evidence that such a subspace is associated with nearly the whole spectral weight generated by 
applying one- and two-electron operators onto m = and x > ground states is provided. It is based on the relative 
amount of such a weight generated by application onto these states of c and av fermion operators. The ground-state 
configurations of such objects were studied in RefJ^. Although part of our results are argued on phenomenological 
grounds, they emerge naturally from the scenario provided by the interplay of symmetry and the suitable quantum- 
object operator description used in this paper. It is found that for the Hubbard model on the square lattice in the one- 
and two-electron subspace only the c and si fermions play an active role. Hence one may neglect the remaining av 
fermion branches considered in Ref.^. The states of such a subspace may have none or one spin- neutral four-spinon 
s2 fermion but such an object has vanishing momentum and vanishing energy. In Ref. 19 important physical quantities 
of the square- lattice quantum liquid introduced here are expressed in terms of c and si fermion energy dispersions 
and velocities. It is confirmed in Sections II-A and II-B that for simple one- and two-electron operators O the leading 
operator term O on the right-hand side of Eq. ^ generates nearly the whole spectral weight. For such operators the 
terms of the general expression Q containing commutators involving the related operator 5 = 5 are found to generate 
very little spectral weight. Hence one can reach a quite faithful representation of such operators by expressing them 
in terms of the c and si fermion operators. 

In this paper strong evidence is found that provided that in the thermodynamic limit N% — > 00 a long-range 
antiferromagnetic order occurs in the ground state of the related isotropic spin- 1/2 Heisenberg model on the square 
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lattice, a similar long range order sets in in that limit in the ground state of the half-filled Hubbard model on the square 
lattice for U/At > 0. However, it is not among our goals providing a mathematical proof that in the thermodynamic 
limit N% — > oo, x — 0, m — 0, and vanishing temperature T = such a long-range antiferromagnetic order sets 
in. Although there is no such a proof, there is a large consensus that it should be so in both the isotropic spin- 1/2 
Heisenberg model on the square lattice^— and in the half-filled Hubbard on the square lattic o n i 26 ~ — . One of our 
aims is though providing useful physical information on how the occurrence of a long-range antiferromagnetic order 
and a short-range spin order at x — and x > 0, respectively, in the thermodynamic limit N% — > oo, m = 0, and 
vanishing temperature T = is related to the properties of the spin effective lattice. Indeed, one of our goals is 
to show that within the operator description used in this paper the study of the effects of hole doping on the spin 
subsystem simplifies. That simplification follows from the independence that the state representations of the spin 
SU(2) symmetry, 77-spin SU(2) symmetry, and hidden U(l) symmetry recently found in Ref. 4 , respectively, have 
within the present description. Such state representations correspond to independent occupancy configurations of the 
spin-1/2 spinons in the spin effective lattice, ry-spin-1/2 7y-spinons in the 77-spin effective lattice, and c fermions in the 
c effective lattice, respectively. Specifically, at m = spin density one of the main effects of hole doping is on the 
number M a = 2S C of sites of the spin effective lattice, which at x = and x > hole concentrations equals and differs 
from that of the original lattice, respectively. (We recall that for x > the spin effective lattice is well defined only 
for N% 1 very large or infinite^.) 

The microscopic processes corresponding to the effective transfer integrals t' = t'{U/At) and t" = t"(U/At) of 
the Hamiltonian (5) expressed in terms of creation and annihilation rotated-electron operators are needed to char- 
acterize the type of order of the square-lattice quantum liquid. This holds for instance concerning the short-range 
incomensurate-spiral spin order considered in Section III-C. The qualitative changes occurring in the spin effective 
lattice due to hole doping are behind processes that destroy long-range antiferromagnetic order not being active and 
being active at x = and for x > 0, respectively. In this paper strong evidence is provided that for < x <C 1 
the ground state has a short-range incomensurate-spiral spin order. The related investigations of Ref. 1 ^ provide ev- 
idence that a spin short-range order exists for < x < x*, whereas for x > the ground state is a spin disordered 
state without short-range spin order. Here x* is a critical hole concentration whose magnitude is for approximately 
U/At > 1.3 an increasing function of U /At. For the intermediate U /At values of interest for the square-lattice quantum 
liquid studies of that reference it reads = 0.23 for U/At ~ 1.3 and a;* = 0.28 for U/At « 1.6. 

Another goal of this paper is defining the symmetry of the m = and x = ground state after its symmetry is 
spontaneously broken in the limit N% — > 00, upon the emergence of the long-range antiferromagnetic order. It is found 
that at x = 0, m = 0, and vanishing temperature T = that state symmetry is broken from 15*0(3) x SO(3) x U(l) for 
Nl > 1 large but finite to [U{2) x U(l)]/Z$ = [SO (3) x 17(1) x U(l)]/Z 2 for N% -> 00. Finally, the U/At dependence 
of the energy order parameters of the x = antiferromagnetic order and < x <C 1 short-range incomensurate-spiral 
spin order are issues also addressed in this paper. 

The paper is organized as follows. The introduction and basic information on the general operator description 
introduced in RefJ^ is given in Section I. In Section II a suitable one- and two-electron subspace is considered. The 
form that the spin and si effective lattices have in it as well as the corresponding elementary excitations are issues also 
addressed in that section. In Section III it is shown that for N% — > 00 our results are consistent with a Mott-Hubbard 
insulating ground state with long-range antiferromagnetic order at half filling and a ground state with short-range 
spin order for a well-defined range of finite hole concentrations. Strong evidence is given that for < x -C 1 and 
intermediate and large values of U/At the short-range spin order has an incommensurate-spiral character. The U/At 
dependence of the energy order parameters is also addressed. Finally, Section III contains the concluding remarks. 

II. THE SQUARE-LATTICE QUANTUM LIQUID: A TWO-COMPONENT FLUID OF CHARGE c 
FERMIONS AND SPIN-NEUTRAL TWO-SPINON si FERMIONS 

In this section we consider a suitable one- and two-electron subspace and study the form that the spin and si 
effective lattices have in it. To achieve our goals, in the following we address the problem of the generation of the 
one- and two-electron spectral weight in terms of processes of c fermions and spinons. The picture that emerges is 
that of a two-component quantum liquid of charge c fermions and spin neutral two-spinon si fermions. It refers to 
the square-lattice quantum liquid introduced in this paper. 



8 



A. The one- and two-electron subspace 

Let |^gs) be the exact ground state for x > and m = and O denote an one- or two-electron operator. Then 
the state, 

d|* G s> = £c,|^>; Cj = (*j\d\9 a s) , (11) 

j 

generated by application of O onto that ground state is contained in the general one- and two-electron subspace. 
This is the subspace spanned by the set of energy eigenstates j) such that the corresponding coefficients Cj are 
not vanishing. Such a subspace must contain all excitations 0\^gs) with the operator O being any of a well defined 
set of operators. It includes the creation and annihilation electron operator and the whole set of simple two-electron 
operators. For x > and m — the c and si fermion occupancies of the ground states |^gs) are found in Rcf. 12 . 

Our goal here is finding what occupancy configurations of the objects of the description of such a reference generate 
a set of excited energy eigenstates {l^)} such that J^j | Cj | 2 ~ 1 for O being a creation and annihilation electron 
operator and all simple two-electron operators. That set of states must contain nearly the whole spectral weight of the 
above one- and two-electron excitations. Their general form is provided below in Section II-E. It refers to a particular 
case of the momentum eigenstates considered in Reff^ 2 -, which in general are not energy eigenstates. Fortunately, the 
set of such states that span the one- and two-electron subspace are found to be energy eigenstates. For an initial x > 
(and x = 0) and m — ground state \^gs) evidence is provided in the following that such states have excitation 
energy lu < 2/i (and to < The inequality w < 2/i applies to some range of finite hole concentrations x > and 
spin density m — 0. In turn, for the fx = and m = absolute ground state with the chemical-potential zero level at 
the middle of the Mott-Hubbard gap the smallest energy required for creation of either one rotated-electron doubly 
occupied site or one rotated-hole doubly occupied site is instead given by the energy scale This justifies why for 
such initial ground state the above inequality uj < 2/z is replaced by w < fj, . 

From the use of the invariance under the electron - rotated-electron unitary transformation of the independent ±1/2 
spinonsi 2 -, one finds that the number L s = [L s __ 1 / 2 + L s ,+i/2] — 25 s of Eq. (|TU1) for a = s of such objects generated 
by application of ./V-electron operators onto a ground state is exactly restricted to the following range, 

L s = [£,, -i/2 + L a , +1/3] - 25 s = 0, 1, 2, TV. (12) 

It follows that for the one- and two-electron subspace only the values L s = [L S} _i/ 2 + L s .+1/2] = 2S S = 0, 1, 2 are 
allowed. Such a restriction is exact for both the model on the square and ID lattices, as well as for any other bipartite 
lattice. 

For a finite number v > 2 of spinon pairs the sv fermions created onto a x > and m = ground state have 
vanishing energy and momentum^. A vanishing spin density m = refers to a vanishing magnetic field H = 0. 
Hence such objects obey the criterion e s „ = 2vhb H = of Eq. (|A10[) of Appendix A, so that they are invariant 
under the electron - rotated-electron unitary transformation. Therefore, for U/4t > they correspond to the same 
occupancy configurations in terms of both rotated electrons and electrons. That reveals that such sv fermions describe 
the spin degrees of freedom of a number 2v = 4, 6, 2C S of electrons. It follows that nearly the whole spectral weight 
generated by application onto a x > and m = ground state of A/"-electron operators refers to a subspace spanned 
by energy eigenstates with numbers in the following range, 

c s c s 

[L s , -1/2 + L s , +1/2 + 2C S - 25 s ] = 0,1,2, ...,AA; C S = ^>A S „; B s = g N sv . (13) 

!/=l !/=l 

Note that owing to the above invariance of the sv fermions with v > 2 spinon pairs, provided that J\f/N% — > and 
B s /N% — s- for N% —> 00 the number B s = N S v is a good quantum number. (This is a generalization of the 
subspace (A) defined in Ref. 12 .) Consistently, the x > (and x = and /j, = 0) and m = ground state and the set 
of excited states of energy to < 2/i (and to < fi°) that span the one- and two-electron subspace considered here have 
no —1/2 77-spinons, rjv fermions, and sv' fermions with v' > 3 spinon pairs so that N vl/ = and N sl/ > — for v' > 3. 
Summation over the set of states that span such a subspace gives indeed Y^j \Cj\ 2 ~ 1 for the coefficients of the one- 

or two-electron excitation 0\^gs) of Eq. (jlip . This holds for O being the electron creation or annihilation operator 
or any of the simple two-electron operators. Consistently, there is both for the model on the ID and square lattices 
an extremely small weight corresponding mostly to states with iV S 3 = 1, which is neglected within the use of the one- 
and two-electron subspace considered here. (Note that while the number restrictions of Eq. (fl2|) are exact, those of 
Eq. (fl^j) are a very good approximation. This is why J2j IQ| 2 ~ 1 rather than £V \Cj\ 2 = 1 for the j summation 
running over the set of N = 1,2 states that span the one- and two-electron subspace as defined here.) 
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TABLE I: The deviations 6N£ = -28S C and numbers 7Vi\ = [2S S + 2N s2 ] of Eq. (H} for the fourteen classes of one- and 
two-electron excited states of the x > and m — ground state that span the one- and two-electron subspace considered in this 
paper, corresponding electron number deviations SNf and 6N±, and independent-spinon numbers L s ,+i/2 and L Sj _i/2 and s2 
fermion numbers N S 2 restricted to the values provided in Eq. flD - The spin S 3 and deviations SS C , SN s i = [5S C — S s — 2iV»2], 
and SNa sl = [8S C + S s ] of each excitation are also provided. 



Thus, according to Eq. (IT3l) the numbers of independent ±1/2 spinons and that of s2 fermions of the excited states 
that span such a subspace are restricted to the following ranges, 

£.,±i/2 = 0,1; iV s2 = 0, for N = 1 , 
2S s + 2N s2 = [i 5 ,_i /2 + £ s ,+i/2 + 2iV s2 ]=0,l,2, for M=2. (14) 

Here N — 1,2 refers to any of the ./V-electron operators O whose application onto the ground state I^gs) generates 
the above excited states, as given in Eq. (jllj) . Furthermore, the number of c fermions and the number of si fermions 
read N c = N = (1 - x) N% and N sl = [N/2 - 2N s2 - S s ] = (1 - x) N%/2 - [2N s2 + S s ], respectively. If in addition 
we restrict our considerations to the LWS-subspace of the one- and two-electron subspace^, then L s _\i 2 — in Eq. 
(Tnj, whereas the values L s +1 / 2 = 0, 1 for N s2 — and J\f — 1 remain valid and in [2S S + 2N s2 ] = 0, 1, 2 one has that 
2.V, /,.. . , ■ for A' 2. 

As shown in Rcf.— , the numbers of independent 77-spinons (a = rf) and independent spinons {a — s) L a ±i/ 2 , the 
total numbers of 77-spinons (a = 77) and spinons (a = s) M a ±\i 2 — [L a< ±\/ 2 + C a ], the number of sites of the spin 
effective lattice iV^ = 2S C) the number of sites of the 77-spin effective lattice iV^ = [N% — 2,5 C ], the number of c 
fermions iV c = 2S C , and the number of c fermion holes = [N% — 2S C ] are good quantum numbers of the Hubbard 
model on the square lattice. The good news is that in the one- and two-electron subspace considered here the numbers 
N% , N s i, Njji, and N s2 are also good quantum numbers of such a model. The number of sites, unoccupied sites, and 
occupied sites of the c and si effective lattices equal those of discrete momentum values, unfilled momentum values, 
and filled momentum values of the c and si bands, respectively. From straightforward manipulations of Eqs. (|A4|) 
and (| A5|) of Appendix A for av = si we find that the number N% = [JV a i + N^] of si effective lattice sites and thus 
of si band discrete momentum values is in general given by, 

c s 

Nl i =[S c + S s + ^-2)N sl/ ]. (15) 

Hence for the one- and two-electron subspace considered in this paper for which N s[/ = for v > 3 one has that 
N% Bl = [S c + S 8 ] is a good quantum number. Consistently, for spin S s =0 that subspace is a subspace (A) as defined 
in Ref^. It follows that N s i = [S c - 2N s2 ], N^ = 2N s2 , and N s2 and hence N% si = [N sl + N^] = S c are good 
quantum numbers. Furthermore, for the remaining spin values S s = 1/2 and S s = 1 such a subspace is a subspace 
(B) as defined in that reference. Hence N si — [S c — S s ], N^ = 2S S , and N% sl — [S c + S s ] are good quantum numbers. 
It then follows from the general properties of the operator description of Ref^ that for the Hubbard model on the 
square lattice in the one- and two-electron subspace the si fermion band microscopic momenta q are good quantum 
numbers as well. 

Moreover, use of the general expressions of si band discrete momentum values and number NjJi — [N% — N s i] 
of si band unfilled momentum values with N% si given in Eq. (|15[) together with the restrictions in the values of 
the numbers of Eqs. (fl"2"j) , (TT5|) . and (fi"4"|) and the exact result proved in Ref, 1 - that one-electron (and two-electron) 
excitations have no overlap with excited states with none and two (and one) si band holes (and hole) reveals that 
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nearly the whole one- and two-electron spectral weight is contained in the subspace spanned by states whose deviation 
SNj? in the number of c band holes and number of si band holes are given by, 

8N£ = -25S C = -SN = 0, =Fl, =f2 , 
N' s \ = S c + S s - N sl = 2S S + 2N s2 

= ±(5N t - 8NJ + 2L St Tl/2 + 2N s2 = L s , _ 1/2 + L s , +1/2 + 2N s2 -0,1,2. (16) 

Here 8N is the deviation relative to the initial ground state in the number of electrons, SNf and 8Ni are those in 
the number of spin-projection f and J, electrons, respectively, iV S 2 is the number of the excited-state s2 fermions, and 
L s .±i/2 is that of independent spinons of spin projection ±1/2. We emphasize that for an initial m — and S s = 
ground state the numbers N s2 = L S] ±1/2 = N' s \ = vanish^. 

The deviations 8N£ and numbers NjJi of Eq. (fTB]) for the fourteen classes of one- and two-electron excited states 
of the x > and m = ground state that span the one- and two-electron subspace, corresponding electron number 
deviations SN-f and 5N±, and independent-spinon numbers L s ^ +1/2 and L Sj —1/2 ancl s 2 fermion numbers N s2 restricted 
to the values provided in Eq. (fM)) are given in Table Q] For N s2 = 1 spin-siglet excited states the s2 effective lattice 
has a single site and the corresponding s2 band a single vanishing discrete momentum value, q — 0, occupied by the s2 
fermion-^. We recall that such a s2 fermion is invariant under the electron - rotated-electron unitary transformation 
and thus has vanishing energy, consistently with the invariance condition of Eq. (jAlQI) of Appendix A for av = s2 
and vanishing magnetic field H = 0. 

As mentioned above, the initial x > and m = ground states of the one- and two-electron subspace have zero 
holes in the si band so that SNj^ = for the excited states. This follows from all M s — 2S C = N spinons being 
confined within the two-spinon bonds of the N s i = M s /2 si fermions. The one- and two-electron subspace is spanned 
by the states of Table fl] generated by creation or annihilation of \8N^\ = 0, 1, 2 holes in the c momentum band and 
N^i = 0, 1, 2 holes in the si band plus small momentum and low energy particle- hole processes in the c band. The 
charge excitations of x > and m — initial ground states consist of a single particle-hole process in the c band 
of arbitrary momentum and energy compatible with its momentum and energy bandwidths, plus small-momentum 
and low-energy c fermion particle-hole processes. Such charge excitations correspond to state representations of the 
global U(l) symmetry and refer to the type of states denoted by "charge" in the table. The one-electron spin-doublet 
excitations correspond to the four types of states denoted by "±l<rel." in Table U where +1 and —1 denotes creation 
and annihilation, respectively, and a =j, 4.. The spin-singlet and spin-triplet excitations refer to the four types of 
states denoted by "singl.spin" and "tripl.spin" in the table. The two-electron excitations whose electrons are in a 
spin-singlet configuration and those whose two created or annihilated electrons are in a spin-triplet configuration 
correspond to the five types of states "±2 fiel." and "±2erel." of that table where +2 and —2 denotes creation and 
annihilation, respectively, of two electrons. 

For the Hubbard model on the square lattice such fourteen types of states are energy eigenstates and nearly exhaust 
the whole one- and two-electron spectral weight. Excited states of classes other than those of the table contain nearly 
no one- and two-electron spectral weight. Such a weight analysis applies to the ID Hubbard model as well. For the 
corresponding quantum liquid describing the Hamiltonian ([1} in the one- and two-electron subspace, the numbers 
2S C , 2S rn 2S S , and — 25f 3 associated with the global 50(3) x 5*0(3) x U(l) symmetry are given by, 

2S c = (l-x)N*; 2S V = -2S* 3 =xN 2 a - 2S S = (l-x) N 2 a -2[N sl +2N s2 ] ; -2Sf - m N 2 a = 25 s -2L s ,_ 1/2 . (17) 

For such a quantum liquid the number of sites of the spin effective lattice = 2S C and corresponding spacing a s 
given in Eq. (|A8[) of Appendix A read, 

Nl = (l-x)N%; a , = ^=^, (1 - x) > l/N 2 a , (18) 

respectively. Its sites refer to those singly occupied by rotated electrons in the original lattice. In turn, the sites of 
the 77-spin lattice also introduced in Ref.— refer to the sites doubly occupied and unoccupied by rotated electrons in 
the original lattice. 

For the Hubbard model in the one- and two-electron subspace the concept of a 77-spin lattice considered in RefJ^ 
is useless. Indeed, for that subspace such a lattice either is empty (x > 0) or does not exist (x = 0). This is because 
the 77-spin degrees of freedom of the states that span that subspace are the same as those of the C v = Ylu=i v ^n v ~ 
(N 2 /2 - S c - Sri) = vacuum \Q V ;N%) of Eq. (jA"3"j) of Appendix A. For states for which S c = N 2 /2 and thus 
— S n ~ the 77-spin lattice does not exist and thus the spin effective lattice is identical to the original lattice. 
This is argued below in Section III-C to be a necessary condition for a spontaneous symmetry breaking and emergence 
of a long-range antiferromagnetic order to occur in the ground state as N 2 — > 00. In turn, for S c < N 2 /2 the 77-spin 
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degrees of freedom correspond to a single occupancy configuration of the N% independent +1/2 ^-spinons. Such 
objects are invariant under the electron - rotated-electron transformation and thus play the role of unoccupied sites 
of the ?7-spin lattice^. Indeed, for the one- and two-electron subspace the corresponding number of 77-spin lattice 
occupied sites C v = Ylu=l v N n v vanishes. Only for states and subspaces for which N 2 ^/N 2 = [1 — 2S C /N%] is finite 
and the inequality < C v < N 2 holds is the concept of a r;-spin effective lattice useful. For a given state, C v + C s is 
the number of sites of the original lattice whose electron occupancy configurations are not invariant under the electron 
- rotated-electron unitary transformation. For the states that span the one- and two-electron subspace the numbers 
C v and C s are given by C n = and C s — [N s i + 2N s2 ], respectively. For x — (and x = 1) one finds N% — and 
N% s = N 2 (and N% — N 2 and N 2 ^ =0), so that there is no ry-spin (and no spin) effective lattice and the spin (and 
77-spin) effective lattice equals the original lattice. 

As further confirmed below, for the one- and two-electron subspace considered here only the c and si fermions play 
an active role. Straightforward manipulations of the general expressions given in Eq. (fT5)) and Eqs. (| A4|) - (| A9|) of 
Appendix A and related expressions provided in Refp^ reveal that for that subspace the number N 2 si of si band 
discrete momentum values, N s i of si fermions, and of si fermion holes are given by, 

Nl, = i N >i + N si\ = [Sc. + S a ] ; N s i = [So -S.- 2N s2 ] ; JV£ = [2S S + 2N s2 ] = 0, 1, 2 , (19) 

respectively. In turn, the corresponding c effective lattice, c momentum band, and c fermion numbers read, 

Nl = [N c + N£] = N 2 ; N c = 2S c =(l-x)N 2 ; N? = xN 2 . (20) 

B. Confirmation for ID that most one-electron spectral weight is generated by processes obeying the 

ranges of Eqs. (jT3|). (flijl . and (fl6|l 

Above the transformation laws under the electron - rotated-electron transformation of the av fermions were used 
to show that nearly the whole spectral weight of the one- and two-electron excitations of the Hubbard model on the 
square lattice is generated by processes obeying the ranges of Eqs. (fl"3|) . (fl4|) . and ([T6]) . Similar results apply to the 
ID model. 

The terms of one- or two-electron operators O on the right-hand side of the first equation of (U]) that generate 
the excitations O^gs) of Eq. (jlip may be expressed in terms of the of c fermion operators, spinon operators, and 
?7-spinon operators given in Eqs. d6j)- (JSJ) - This is done by use of the operator relations provided in Eq. (|A1|) of 
Appendix A. Concerning the contributions to the O expression provided in Eq. (|4]) containing commutators involving 
the operator S = —(t/U) [T+i — T_i] + 0(t 2 /U 2 ), to fulfill such a task one takes into account that independently of 
their form, the additional operator terms 0(t 2 /U 2 ) of higher order are products of the kinetic operators To, T+i, and 
Tli of Eq. ©. 

From such an analysis one finds that the elementary processes associated with the one- and two-electron subspace 
number value ranges of Eqs. (1131) . (|14[) . and (JT5J) are fully generated by the leading-order operator O. In turn, the 
processes generated by the operator terms containing commutators involving the operator S refer to excitations whose 
number value ranges are different from those of Eqs. (fT3|) . ([M)) . and (|T6|). This confirms that such processes generate 
very little one- and two-electron spectral weight, consistently with the exact number restrictions of Eq. (|12j) . the 
approximate number restrictions of Eq. (|13l) , and the results of Ref J^. 

For the Hubbard model on the ID lattice also often considered in the studies of that reference, the spectral-weight 
distributions can be calculated explicitly by the pseudofermion dynamical theory associated with the model exact 
solution, exact diagonalization of small chains, and other methods. From the use of the same arguments as for 
the model on the square lattice one finds that at ID the one- and two-electron subspace considered in this paper 
corresponds to the same number deviations and numbers as for the square lattice. In addition, the relative one- 
electron spectral weight generated by different types of microscopic processes can be studied by means of the above 
methods. That program is fulfilled in Ref. 32 . The results of that reference confirm the dominance of the processes 
associated with the number value ranges provided in Eqs. and (|16[) . They refer spacifically to operators O and 
O that are electron and rotated-electron, respectively, creation or annihilation operators. Such studies confirm that 
the operator O generates all processes associated with the number value ranges of Eqs. (|14|) and ([T6|) and number 
values of Table HI In addition, it also generates some of the non-dominant processes. That is confirmed by the weights 
given in Table 1 of Ref. 32 , which correspond to the dominant processes associated with only these ranges. The small 
missing weight refers to excitations whose number value ranges are not those of Eqs. (|14j) and (|16j) but whose weight 
is also generated by the operator O. Indeed, that table refers to U /At — > 00 so that O = O and the operator terms of 
the O expression provided in Eq. (j4} containing commutators involving the operator S vanish. 
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On the other hand, for finite values of U/At all dominant processes associated with the number value ranges of Eqs. 
(1141) and (|16|) and number values provided in Table U are generated by the operator O. In turn, the small spectral 
weight associated with excitations whose number value ranges are different from those are generated both by that 
operator and the operator terms of the O expression of Eq. (j4|) containing commutators involving the operator S. 
For the model on the ID lattice the small one-electron spectral weight generated by the non-dominant processes is 
largest at half filling and U /At ~ 1. 

For x > the one- and two-electron subspace is spanned by states with vanishing rotated-electron double occupancy. 
This holds both for the Hubbard model on the ID and square lattices. Generalization of the results to the range x < 
reveals that then such a subspace is spanned by states with vanishing rotated-hole double occupancy. That property 
combined with the particle-hole symmetry explicit at x — and \l = 0, implies that the relative spectral- weight 
contributions from different types of one-electron addition excitations given in Fig. 2 of Ref. 32 for the ID model 
at half filling leads to similar corresponding relative weights for half-filling one-electron removal. Analysis of that 
figure confirms that for the corresponding one-electron removal spectrum the dominant processes associated with the 
number value ranges of Eqs. (fT4"]) and (fTf?)) and number values given in Table U refer to the states called in figure 
1 holon - 1 si hole states. Their minimum relative weight of about 0.95 is reached at U/At « 1. For other hole 
concentrations x > and values of U/At the relative weight of such states is always larger than 0.95, as confirmed 
from analysis of Figs. 1 and 2 and Table 1 of Refi^ 2 -. 

For the Hubbard model on the square lattice the explicit derivation of one- and two-electron spectral weights is a 
more involved problem. The number value ranges of Eqs. (|13[) , (|14[) , and (|16l) and number values provided in Table 
U for the one- and two-electron subspace also apply, implying similar results for the relative spectral weights of one- 
and two-electron excitations. Consistently and as mentioned above, there is an exact selection rule valid both for the 
Hubbard model on a ID and square lattices that confirms that = 1 and = 0, 2 for one-electron excitations 
and two-electron excitations, respectively. It follows from the expression provided in Eq. (|A7[) of Appendix A for the 
quantum number = e i7rN ^ = e i7rN = ±1. It reveals that for SN = ±1 (and SN = ±0 and SN = ±2) excited 
states the number of holes in the si band must be always an odd (and even) integer. Here SN is the deviation 
in the electron number N under a transition from a x > and m = ground state to such excited states. This 
implies that one-electron (and two-electron) excitations do not couple to excited states with two holes (and one hole) 
in the si band. Indeed, for such excitations one has that e mSN = — 1 (and e" SN = +1), so that e l7rSN * 1 — — 1 (and 

e fer«JV* =+1). 



C. The spin and si effective lattices for the one- and two-electron subspace 

According to the restrictions and numbers values of Eqs. (|14p and (fTBT) and Table HI the states that span the one- 
and two-electron subspace may involve none or one s2 fermion. As confirmed in Ref.—, it is convenient to express the 
one- and two-electron excitation spectrum relative to initial x > and m = ground states in terms of the deviations 
in the numbers of c effective lattice and si effective lattice unoccupied sites. Those are given explicitly in Eq. (I16[) and 
TableU Note that the number of si fermions provided in Eq. (|T9|) can be written as N s i — [(1 — x)N 2 /2 — S s — 2N S 2] 
where S s = for N s2 = 1 and S s = 0, 1/2, 1 for N s2 = 0. 

As discussed above, for N s2 = 1 spin-singlet excited energy eigenstates the single s2 fermion has vanishing energy 
and momentum and consistently with Eq. (IA10[) of Appendix A, for vanishing magnetic field if = it is invariant 
under the electron - rotated-electron unitary transformation. Therefore, the only effect of its creation and annihilation 
is in the numbers of occupied and unoccupied sites of the si effective lattice. Its creation can then be accounted 
for merely by small changes in the occupancies of the discrete momentum values of the si band, as discussed below. 
Hence the only composite object whose internal occupancy configurations in the spin effective lattice are important 
for the physics of the Hamiltonian ([I} in the one- and two-electron subspace is the spin- neutral two-spinon si fermion 
and related spin-singlet two-spinon si bond particl e 12 ' 31 . 

It is confirmed below that for the Hubbard model in the one- and two-electron subspace and alike for the s2 fermion, 
the presence of independent spinons is felt through the numbers of occupied and unoccupied sites of the si effective 
lattice. In turn, the number of independent +1/2 ?7-spinons equals that of the unoccupied sites of the c effective 
lattice. Therefore, when acting onto that subspace, the Hubbard model refers to a two-component quantum liquid 
that can be described only in terms of c fermions and si fermions. 

For x > and states belonging to the one- and two-electron subspace the spin effective lattice has a number of sites 
N 2 = (1 — x) N 2 . Its value is smaller than that of the original lattice. The lattice constant a s provided in Eq. (|A8|) 
of Appendix A for a = s reads a s « a/\J\ — x for such states, as given in Eq. (fT5|). It is such that the area L 2 of the 
system is preserved. Any real-space point within the spin effective lattice corresponds to the same real-space point 
in the system original lattice. Except for a suitable phase factor, a local si fermion has the same internal structure 
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as the corresponding si bond-particl o 12 ' 19 . The si fermion spinon occupancy configurations considered in Rcf. 31 are 
expected to be a good approximation provided that the ratio N% /N% and thus the electronic density n = (1 — x) are 
finite in the limit N% — > oo. This is met for the hole concentration range x £ (0, x*) where x* < 1 for which according 
to the studies of Ref^ the maximum magnitude of the si fermion spinon pairing energy is finite. 

Within the N% 3> 1 limit that the description used in the studies of this paper refers to, there is for the states that 
span the one- and two-electron subspace commensurability between the real-space distributions of the ~ N a i 
sites of the si effective lattice and the N% ~ 2N s i sites of the spin effective lattice. For (1 — x) > 1/N% and iV| 3> 1 
the spin effective lattice has N% = (1 — x) N% sites and from the use of the expression given in Eq. (|T5j) for the 
number of si effective lattice sites N% and Eq. (|A9| of Appendix A for the corresponding spacing a s i we find, 



2S 
(l-x)N% 



— 2 9 1 \ 1 

( ' -2(1^) N§)^ a " ^ = °'2' 1 - < 21 > 



For = states such as the x > and m = ground states the square spin effective lattice has two well-defined 
sub-lattices, which we call sub-lattice A and B, respectively. As discussed in Refi^i, for =1,2 states the spin 
effective lattice has two bipartite lattices as well, with one or two extra sites corresponding to suitable boundary 
conditions. The two spin effective sub-lattices have spacing a s i = \/2a s . The fundamental translation vectors of the 
sub-lattices A and B read, 

a s i = ^=(e Xl + e X2 ) , b sl = --^=(e Xl - e X2 ) , (22) 

respectively. Here e Xl and e X2 are the unit vectors and x\ and X2 Cartesian coordinates. As confirmed in Ref^I, the 
vectors given in Eq. (f2"2l are the fundamental translation vectors of the si effective lattice. 

In the case of x > 0, m = 0, and = ground states whose si momentum band is full and all N% sites of the 
si effective lattice are occupied we consider that the square root N as of the number N% of sites of the spin effective 
lattice is an integer. It then follows that the spin effective lattice is a square lattice with N as x N aB sites. Thus the 
square root N asl of the number N% of sites of the si effective lattice cannot in general be an integer number yet 
N% is. However, within the N% 3> 1 limit considered here we use the notation N% for the number of sites of the si 
effective lattice. 



D. The quantum liquid of c fermions and si fermions: Why only such objects play an active role? 

For the one- and two-electron subspace considered in this paper the number N% of sites of the si effective lattice 
and si band discrete momentum values, N s i of si fermions, and of si fermion holes have expressions given in Eq. 
(jTSJ). The corresponding numbers of the c effective lattice and c band are provided in Eq. (|2T)]). For that subspace the 
si band is either full and thus is filled by N s i = AT? — 2S C = si fermions or has one or two holes. Furthermore, 
one-electron and two-electron excitations have no overlap with excited states with two holes and one hole in the si 
band, respectively. Specifically, excited states with a single hole in the si band correspond to one-electron excitations 
and those with — 0, 2 holes in that band refer to two-electron excitations, as given in Table |TJ Excited states 
with = 3 (and = 4) holes in the si momentum band correspond to very little one-electron (and two-electron) 
spectral weight and are ignored within the use of the one- and two-electron subspace considered here. 

We now justify why the square-lattice quantum liquid corresponding to the Hubbard model on the square lattice 
in the one- and two-electron subspace may be described only by c and si fermions on their c and si effective lattices, 
respectively. According to the number value ranges of Eqs. (jT4j) and (| 16[) and number values provided in Table HI the 
one- and two-electron subspace is spanned by excited states having either none A^ S 2 = or one N S 2 = 1 spin-neutral 
four-spinon s2 fermion. = 1 spin-singlet excited states have no independent spinons. One then finds from the use 
of Eq. (|A5|) of Appendix A for av — s2 that N^ 2 = so that such states have no holes in the si momentum band and 
thus N% s2 = 1. This means that for such states the s2 fermion occupies a s2 band with a single vanishing momentum 
value. Since the s2 fermion under consideration has both vanishing momentum and energy and is invariant under 
the electron - rotated-electron unitary transformation, the only explicit effect of its creation is onto the numbers of 
occupied and unoccupied sites of the si effective lattice and corresponding numbers of si band si fermions and si 
fermion holes. Specifically, according to the expressions provided in Eq. ([!"§]) and number values of Table QJ the 
deviations SS C — 5S S — and SN S 2 — 1 generated by a state transition involving creation of one s2 fermion lead to 
deviations in the number of si fermions and si fermion holes given by SN s i = —2SN S 2 — —2 and SN^ = 28N S 2 = 2, 
respectively. 
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Moreover, the ranges of Eqs. (|14p and (fT5f and number values of Table U confirm that such iV S 2 = 1 excited states 
have zero spin, S s — 0. According to Eq. (fT9)) , the number of holes in the si band is Nj^ — 2N S 2 — 2 for such states, 
in contrast to N' s \ = for the initial ground state. In turn, the number A 2 si of sites of the si effective lattice remains 
unaltered. Following the annihilation of two si fermions and creation of one s2 fermion, two unoccupied sites appear 
in the si effective lattice. As a result two holes emerge in the si band as well. The emergence of these unoccupied 
sites and holes involves two virtual processes where (i) two si fermions are annihilated and four independent spinons 
are created and (ii) the latter independent spinons are annihilated and the s2 fermion is created. 

Hence the only explicit net effect of the creation of a single vanishing-energy and zero-momentum s2 fermion is 
the annihilation of two si fermions and corresponding emergence of two holes in the si band and two unoccupied 
sites in the si effective lattice. Therefore, in the case of the one- and two-electron subspace one can ignore that 
object in the theory provided that the corresponding changes in the si band and si effective lattice occupancies are 
accounted for. Within neutral si fermion particle-hole processes of transitions between two excited states with a 
single s2 fermion, two of the four spinons of such an object are used in the motion of si fermions around in the si 
effective lattice. Indeed, such two spinons play the role of unoccupied sites of that lattic o 12 ' 31 , consistently with the 
expression N> s \ = 2N s2 given in Eq. flU]) for 2S S = 0. 

Also the L s = 2S S independent spinons play the role of unoccupied sites of the si effective lattice. Again, this is 
consistent with the expression — L s = 2S S = 1,2 provided in Eq. (IT9)) for the number of unoccupied sites of the 
si effective lattice and of si fermion holes of the corresponding N S 2 = excited states. As given in Eqs. (|T4"f and (flB]) 
and Table HI the one- and two-electron subspace L s allowed values are L s — [L s t -x/2 ± L s ,+1/2] = 25 s = 0, 1, 2. For 
25 s = 1,2 one has that N S 2 = 0. Now in contrast to creation of a single s2 fermion, a deviation S2S S = 1,2 generated by 
a transition from the ground state to such 2S S = 1,2 excited states may lead to deviations in the numbers of occupied 
and unoccupied sites of the si effective lattice and corresponding si fermion and si fermion holes that do not obey 
the usual equality SN s i = —SN^. Indeed, in the present case 25S C = ±1 for SN^ = 25S S = 1 and 25S C = 0, ±2 for 
SNg ± = 26 S s = 2. Hence according to the expressions provided in Eq. (fTT)|). such deviations lead to deviations in the 
numbers of occupied and unoccupied sites of the si effective lattice and corresponding numbers of si fermions and si 
fermion holes. Those read 8N s i = [SS C — 8S S ] and SNjji — 52S S , respectively. It follows that the total number of sites 
and thus of discrete momentum values of the si band may change under such transitions. This leads to an additional 
deviation 6 A 2 ^ = [SS c + 8S a ]. As given in Table HI for one-electron excited states one has that SN^i — 28S S = 1 
and 26S C = ±1. As a result, 8N sl = ±1/2 - 1/2 = -1,0 and 6N% si = ±1/2 + 1/2 = 0, -1. In turn, for N s2 = 
two-electron excited states one has SNjfr = 2SS S = 2 and 2SS C = 0, ±2 . Thus SN s i = -1, (±1 - 1) = -2,-1, and 
8Nl tl = 1, (±1 + 1) =0,1, 2. 

For the si fermion operators f~ sl and /<f jS i, excitations that involve changes 8N^ si = [8S C ± 8S S ] in the number 
of sites and discrete momentum values of the si effective lattice and si band, respectively, correspond to transitions 
between different quantum problems. Indeed, such operators act onto subspaces spanned by neutral states, which 
conserve S c , S s , and N% sl - In turn, the generator of a non-neutral excitation is the product of two operators. The 
first operator makes small changes in the si effective lattice or corresponding si momentum band. Such changes 
follows the above deviations <5A 2 si = [8S C ± 8S S ]. The second operator is a si fermion operator or a product of such 
operators appropriate to the excited-state subspace. 

Also the vanishing momentum and energy L v+ i/2 = x A 2 independent ±1/2 ?7-spinons are invariant under the 
electron - rotated-electron unitary transformation. Their creation or annihilation may be accounted for by small 
suitable changes in occupancies of the c effective lattice and c momentum band. For x > and the one- and two- 
electron subspace considered here such independent ±1/2 77-spinons correspond to a single occupancy configuration 
associated with the ry-spin vacuum |0 ?) ; A 2 ^) of Eq. (|A3I) of Appendix A. In turn, the degrees of freedom of the 

rotated-electron occupancies of such 2; A 2 sites of the original lattice associated with the U(l) symmetry refer to 
the unoccupied sites of the c effective lattice of Eq. ([20]) and corresponding c band holes. Hence the number 
^77, +1/2 = x A 2 of independent ±1/2 ry-spinons equals that 2V£ = a; A 2 of unoccupied sites of the c effective lattice 
and corresponding c band holes. This confirms that the deviations 8L v +1 /2 = (&e)A 2 originated by creation and 
annihilation of independent ±1/2 ?7-spinons within the one- and two-electron subspace has no effects on the physics 
other than the corresponding deviation SN^ = (8x) iV 2 in the number of unoccupied sites of the c effective lattice and 
c band holes. 

Spin-singlet excitations generated by application onto a m — and x > initial ground state of the operator 
/o s2 fq,si fq 'si where q and q' are the momenta of the two emerging si fermion holes are neutral states which conserve 
S c , S s , and A 2 si . (See Table HI) The implicit role of the s2 fermion creation operator /J s2 is exactly canceling the 
contributions of the annihilation of the two si fermions of momenta q and q' to the commutator [q s i Xl -, Qsix 2 ] °f the 
si translation generators in the presence of the fictitious magnetic field B s i of Eq. (IA2|) of Appendix A. This ensures 
that the overall excitation is neutral. Since the s2 fermion has vanishing energy and momentum and the si band and 
its number A 2 si of discrete momentum values remain unaltered, one can effectively consider that the generator of 
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such an excitation is fg- s i fq> s \ and omit the s2 fermion creation operator. Its only role is ensuring that the overall 
excitation is neutral and the two components of the si fermion microscopic momenta can be specified. It follows 
that for the one- and two-electron subspace the operators fe al /^ >s i, f\i sX /^ )Sl , f\, s i fq',si> an d fa s i f\> a \ generate 
neutral excitations. 

In summary, when acting onto the one- and two-electron subspace considered in Section I-A, the Hubbard model on a 
square lattice refers to a two-component quantum liquid described in terms of two types of objects on the corresponding 
effective lattices and momentum bands: The charge c fermions and spin-neutral two-spinon si fermions. The one- and 
two-electron subspace can be divided into smaller subspaces that conserve S c and S s . Those are spanned by states 
of general form given below in Section II-E. When expressed in terms of c and si fermion operators, the Hubbard 
model on a square lattice in the one- and two-electron subspace is the square-lattice quantum liquid further studied 
in ReffiS. The presence of independent +1/2 spinons or of a composite s2 fermion is accounted for by the values of 
the occupied and unoccupied sites numbers of the si effective lattice and corresponding si fermion and si fermion 
holes. In turn, the number of independent +1/2 ?7-spinons equals that of the unoccupied sites of the c effective lattice 
and c band holes. Otherwise, the presence of vanishing momentum and energy independent spinons or of a single 
spin- neutral four-spinon s2 fermion as well as that of independent +1/2 ?7-spinons has no explicit direct effects on 
the physics. This property follows from all such objects being invariant under the electron - rotated-electron unitary 
transformation^. 

The quantum-liquid c fermions are ?7-spinless and spinless objects without internal degrees of freedom and structure 
whose effective lattice is identical to the original lattice. For the complete set of U /4i > energy eigenstates that 
span the Hilbert space the occupied sites (and unoccupied sites) of the c effective lattice correspond to those singly 
occupied (and doubly occupied and unoccupied) by the rotated electrons. The corresponding c band has the same 
shape and momentum area as the first Brillouin zone. 

In contrast, the quantum-liquid composite spin-neutral two-spinon si fermions have internal structure and the 
definition of the si effective lattice in terms of both the original lattice and the spin effective lattice as well as the 
spinon occupancy configurations that describe such objects is for the one- and two-electron subspace a more complex 
problem 31 . It is simplified by the property of the states that span such a subspace that the corresponding si effective 
lattice has none, one, or two unoccupied sites. 



E. The c and si fermion momentum values and the energy eigenstates 

Here we provide the specific form that the momentum energy eigenstates considered in Refj^ 2 - have in the one- and 
two-electron subspace. Such states refer to a complete set of states in the full Hilbert space. In general they are not 
energy eigenstates. Fortunately, in the one- and two-electron subspace such momentum energy eigenstates are as well 
energy eigenstates. This confirms the usefulness of the square-lattice quantum liquid that refers to that subspace. 

The si band discrete momentum values qj where j = 1, N% are the conjugate of the real-space coordinates fj 
of the si effective lattice for which also j = 1, N% . The same applies to the c band discrete momentum values qj 
and the c effective lattice real-space coordinates fj where in both cases j = 1, N%. (The latter lattice is identical to 
the original lattice.) The c translation generators q c commute with both the Hamiltonian and momentum operator 
for the whole Hilbert spaced. This is why the c band discrete momentum values are good quantum numbers. In 
turn, the si translation generators q sl in the presence of the fictitious magnetic field B s \ of Eq. (IA2[) of Appendix 
A do not commute in general with the Hamiltonian of the Hubbard model on the square lattice. Combining the 
results of Rcf. 1 — with the specific properties of that model in the one- and two-electron subspace reveals that in the 
neutral subspaces of such a subspace the si translation generators q sl do commute with both the Hamiltonian and 
momentum operator. This is why for the present square-lattice quantum liquid the si fermion discrete momentum 
values q = [(fei, c^] are good quantum numbers and thus are conserved. The c and si translation generators reacts, 

<f c = QN C (q) ; £ x = £ q N sl (q) . (23) 
<? <? 

Here N c (q) and N s ±(q) are the momentum distribution-function operators, 

NM = fi c fo c ; N.i(& = fl al fo.u (24) 

respectively. For the Hubbard model on the square lattice in the one- and two-electron subspace the expression of the 
momentum operator simplifies. It reads, 



P = q c + q sX 



(25) 
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Indeed, we recall that in it the c2 fcrmion, independent ±1/2 spinons, and independent +1/2 77-spinons have vanishing 
momentum. 

Since in contrast to the c fcrmions, the si fermions have internal structure, how is the si fermion momentum q 
related to the two underlying spinons? Independent spinons carry no momentum and are invariant under the electron 
- rotated-electron unitary transformation 12 . On the other hand, within the LWS representation of the spin SU(2) 
algebra^, the spin-down spinon of the spin-singlet two-spinon the si fermion of momentum q carries momentum q 
and its spin-up spinon carries momentum — q. In turn, within the highest-weight state (HWS) representation of that 
algebra, its spin-down spinon carries momentum — q and its spin-up spinon carries momentum q. Alike in Refi^ 2 -, 
here we use the LWS representation, so that the spin-singlet two-spinon si fcrmions of momenta (i) q and (ii) — q 
involve (i) a spin-down spinon of momentum q and a spin- up spinon of momentum — q and (ii) a spin-down spinon of 
momentum — q and a spin-up spinon of momentum q, respectively. 

Within the LWS representation, a one-electron removal excitation breaks the spin-singlet spinon pair of a si fermion 
before the annihilation of a spin-down electron. The spin-down spinon of momentum q is then removed within the 
electron. The uncompensated spin- up spinon momentum — <f is associated with that of a hole emerging in the si band 
at momentum q. Indeed, the latter spinon decays into that momentum — q si band hole and a vanishing-momentum 
spin-up independent spinon. Hence one-electron excitations break spinon bond pairs whose spinons had momenta q 
and — q corresponding to their relative motion in the pair. 

In turn, spinon pair breaking under spin excitations or excitations involving removal of two electrons with the same 
spin projection may introduce an extra momentum contribution that corresponds to the motion of the center of mass 
of the broken spin-singlet spinon pair. Under such excitations there emerge two holes in the si band at momenta q' 
and q" of the general form, 

q+ Uq. (26) 

q' + q" . (27) 

where q corresponds to the spinon relative motion in the pair and Sq refers to the motion of the center of mass of the 
spinon pair. 

That for the square-lattice quantum liquid both the c and si fermion discrete momentum values are good quantum 
numbers confirms the suitability of the present description in terms of occupancy configurations of the c and si 
effective lattices. Indeed, the c and si band discrete momentum values are the conjugate of the real-space coordinates 
of the c and si effective lattice, respectively. Are the approximations used in the construction of the si effective lattice 
inconsistent with the si band discrete momentum values being good quantum numbers? The answer is no. Indeed, 
such approximations concern the relative positions of the j — 1, ...,N% sites of the si effective lattic o 12 ' 19 ! 31 . They 
are only directly related to the shape of the si band boundary. They do not affect the si band discrete momentum 
values being good quantum numbers. At x = the spin effective lattice is identical to the original square lattice 
and the si effective lattice is one of its two sub-lattices. Consistently, at x = and m = the boundary of the si 
momentum band is accurately known. Indeed, then the si band coincides with an antiferromagnetic reduced Brillouin 
zone of momentum area 27r 2 such that \q Xl \ + \Qx 2 \ —• n ~- I* 1 turn, it is known that for x > and m = the the si 
band boundary encloses a smaller momentum area (1 — x)27r 2 yet its precise shape remains an open issue. The related 
problems of the c and si momentum bands and corresponding energy dispersions and velocities are studied in Ref^ 9 -. 

For a number of sites N% S> 1 very large but finite that the description used in the studies of this paper refers to a 
m = ground state is both for x = and x > a spin-singlet stated 2 -. For m — and x = this agrees with a theorem 
introduced and proved in Ref£. The corresponding one- and two-electron subspace considered in this paper may be 
divided into a well-defined set of smaller subspaces spanned by neutral states. Such states conserve the eigenvalue 
S c = [1/2] (1 — x)N% of the generator of the hidden U(l) symmetry and the spin S s = 0, |, 1 and thus conserve as well 
the number of sites of the si effective lattice N% sl = [S c + S s ]. The set of energy eigenstates that span such subspaces 
are particular cases of the general momentum eigenstates \&u/4t) studied in Ref.— . The use of the general expression 
of such states leads to the following corresponding general form for the energy eigenstates \^u/4t) = \®u/4t) that span 
the one- and two-electron subspace considered here, 

rst) L «.-v« Vr 172 

\*V/4t) - /7T- \*LWS-,v/4t); C s = 5 Ls _ 1/2 , „ + [] I [L s + 1 - / ] = 1, 2, 4 . (28) 

vcs 1=1 

The LWS appearing in this equation reads, 

l*ws;^} = [|0,;A3^ 4^1 = ^4,^ flc = ^4, e V. (29) 

q' q 



Here, 



q' = q + ^q; q" 
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Here .Fl, s1 and J-1 c are the creation operators of aU/At — > oo si fermion of momentum q' and c fermion of momentum 
q, respectively^. Moreover, \0 V ; N 2 ^) is the 77-spin SU(2) vacuum associated with N 2 independent +1/2 ?7-spinons, 

|0 S ; N 2 J is the spin SU(2) vacuum associated with independent +1/2 spinons, and \GS C ; 0) is the c U (1) vacuum. 
Such three vacua are invariant under the electron - rotated-electron unitary transformation, refer to the model global 
[SU(2) x SU(2) x U(l)]/Z$ = 50(3) x 50(3) x (7(1) symmetry*, and appear in the theory vacuum of Eq. CM) of 
Appendix A. (In that equation, \GS C ;2S C ) = Uq fl c \GS c ;0) .) 

The more general states \$>u/4t) considered in Ref. 12 involve occupancy configurations of the remaining r/i' fermion 
and v > 1 sv fermion branches and —1/2 77-spinon occupancies absent in the expressions given in Eqs. (1281) and (1291) . 
Importantly, the results of that reference concerning the subspaces A and B defined in it confirm that the states of 
form ([28]) and (|29|) are indeed energy eigenstates. Since they span all subspaces of the one- and two-electron subspace 
that conserve S c and S s , they span the latter subspace as well. In contrast, the momentum eigenstates generated 
by simple occupancy configurations of c and av fermions and independent 77-spinons and spinons of the larger set of 
states {| < I ) (7/4t)} considered in Ref.— are not in general energy eigenstates. As justified in that reference, the energy 
eigenstates are superpositions \^u/u) = Xz ^2 \^u/4t-i) 01 a set 01 such states {\^u/4t-i}} with the same momentum 
eigenvalue. We recall that states with a single s2 fermion have also the general form provided Eqs. (|2"8")l and (f2T)|). As 
discussed above, the presence of that vanishing-energy, vanishing-momentum, and spin-neutral four-spinon object is 
accounted for the values of the numbers [N^ - 2S S ] = [S c - S s - N sl ] = 2N s2 = 0, 2 of Eq. 1(15]). 

The energy eigenstates \&u/4t) °f general form (|2"5) that span the one- and two-electron subspace have numbers 
N s2 = A 2 s9 = 0, 1 and N s i w N 2 gi such that N' s \ = [N 2 gi - N sl ] = 0, 1,2. Hence according to the results of Ref. 12 

the lattice spacing a s % ~ \[2a s = yj2j(A — x) a of Eq. (|2Tj) is directly related to the fictitious magnetic-field length 
l s i associated with the field of Eq. (|A2I) of Appendix A. Indeed, in that subspace one has that (np. tS i) « 1 and such 

a fictitious magnetic field reads B s \(fj) ss <&o Y]jiy.j S(rji — fj) e X3 . It acting on one si fermion differs from zero only 
at the positions of other si fermions. In the mean field approximation one replaces it by the average field created 
by all si fermions at position rj. This gives, B s i(fj) w $ n s i(fj) e X3 « 4> [N^ 3l /L 2 ] e X3 = [$o/a 2 i] e X3 . One then 
finds that the number iV 2 si of the si band discrete momentum values equals [B sl L 2 ]/$ - In addition, the si effective 

lattice spacing a s i is expressed in terms to the fictitious magnetic-field length l s i « a/ \/ir(l — x) as a 2 sl — 2irl 2 sl . 
This is consistent with each si fermion having a flux tube of one flux quantum on average attached to it. 

As further discussed in RefJ^, for the present one- and two-electron subspace the si fermion problem is then related 
to the Chern-Simons theory^. Indeed the number of flux quanta being one is consistent with the si fermion and si 
bond-particle wave functions obeying Fermi and Bose statistics, respectively. Hence the composite si fermion consists 
of two spinons in a spin-singlet configuration plus an infinitely thin flux tube attached to it. Thus, each si fermion 
appears to carry a fictitious magnetic solenoid with it as it moves around in the si effective lattice. 

III. LONG-RANGE ANTIFERROMAGNETIC ORDER AND SHORT-RANGE 
SPIRAL-INCOMMENSURATE SPIN ORDER FOR x = AND < x < 1, RESPECTIVELY 

Here we profit from the rotated-electron description used in the studies of this paper to address issues related to 
the occurrence for m = 0, zero temperature T = 0, and N 2 -> 00 of a long-range antifcrromagnetic order and a 
short-range spiral- incommensurate spin order at x — and for < x <C 1, respectively. The emergence as A^ 2 — > 00 
of a long-range antiferromagnetic order in the x = and m = ground state is associated with a spontaneous 
symmetry breaking. We argue that a condition for emergence of such a long-range antiferromagnetic order is that 
the spin effective lattice is identical to the original lattice. Such a condition is not met for small hole concentrations 
< x < 1. 

A. Extension of the Mermin and Wagner Theorem to the half-filled Hubbard model for U /At > 

It is well known that for U /At ^> 1 the spin degrees of freedom of the half- filled Hubbard model on a square lattice 
may be described by an isotropic spin-1/2 Heinsenberg model on a square lattice. It follows that the Mermin and 
Wagner Theorem^ is valid for the former model at half filling and U /At ^> 1. The theorem states that then there is 
no long-range antiferromagnetic order for finite temperatures T > and N 2 — >• 00. 

Let us provide evidence that the Mermin and Wagner Theorem applies to the half-filled Hubbard model on a 
square lattice for all values U/At > 0. The possibility of such an extension to U/At > is strongly suggested by 
evidence involving the transformation laws of the spin configurations under the electron - rotated-electron unitary 
transformation. We recall that in terms of the rotated electrons as defined in Section I and RefJ^, the occupancy 
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configurations that generate the energy eigenstates Y&u/it) = V^l^oo) are the same for all finite values U/At > 
0. Moreover, such rotated-electron occupancy configurations equal those that generate the corresponding energy 
eigenstates l^oo) in terms of electrons in the U/At — > oo limit. 

The rotated-electron configurations that generate the energy eigenstates Y&u/At) that span the one- and two-electron 
subspace defined in this paper are much more complex than those associated with the simple form (|28[) and (|29l) in 
terms of c and si fermion operators. Indeed, the expression of their generators is in terms of rotated-electron 
creation operators is an involved problem. This follows from the expression of the spin- neutral two-spinon si fermion 
operators not being simple in terms of the rotated-electron operator a 19 ' 31 . This reveals that the electronic occupancy 
configurations that in the U/At — » oo limit generate such energy eigenstates correspond to an involved problem as 
well. For U /At — ¥ oo the c fermion holes, spinons, and ?7-spinons are the "quasicharges" , spins, and pseudospins, 
respectively, of Reb 3 ^. For the Hubbard model on the square lattice in the one- and two-electron subspace the energy 
bandwidth of the si fermion dispersion vanishes for U/At —¥ oo and the c fermion dispersion has in that limit the 
simple form e c (q) = — 2£^i=i [ cos (9^i) — C0S ( c lFcx i )} in terms of the c band momentum components^ 9 -. Here qpcxi 
where % = 1, 2 are the components of the c Fermi momentum qp c defined in Reb^ 9 -. 

The use below of the following two properties provides strong evidence that the Mermin and Wagner Theorem 
holds for the half-filled Hubbard model on the square lattice for U/At > 0: I) The x = and m = absolute ground 
state is in the limit N% — > oo invariant under the electron - rotated-electron unitary transformatio n 12 ! 19 . Hence the 
occurrence for N% — > oo of long-range antiferromagnetic order as U/At — > oo, associated with that of the isotropic 
spin-1/2 Heisenberg model, implies the occurrence for iV 2 — » oo of that long-range order for U/At > as well; II) 
Since in terms of rotated electrons single and double occupancy are good quantum numbers for U/At > 0, the rotated- 
electron occupancy configurations that generate the energy eigenstates are more ordered than the corresponding 
electron occupancy configurations. It follows that the lack of long-range antiferromagnetic order in terms of the spins 
of the rotated electrons implies as well a lack of such an order in terms of the spins of the electrons whose occupancy 
configurations generate the same states. 

The rotated-electron operator description of RefJ^ has been constructed to inherently the electron occupancy 
configurations that for U/At — > oo generate an energy eigenstate \^oo) being identical to the rotated-electron occupancy 
configurations that for U/At > generate the energy eigenstates \^u/it) = V^l^oo) belonging to the corresponding V 
tower. Hence concerning the original-lattice rotated-electron occupancies, the Mermin and Wagner Theorem applies 
to all finite values of U/At > 0: That for the occupancy configurations of the rotated-electron spins there is no 
long-range antiferromagnetic order for temperatures T > and TV 2 — > oo is an exact result. The above property II 
then implies that the lack of long-range antiferromagnetic order of the rotated-electron spins for U/At > and T > 
implies a similar lack of such an order for the spins of the original electrons. Indeed, for finite values of U/At the spin 
occupancy configurations are more ordered for the rotated electrons than for the electrons. 

This is also consistent with the emergence of a long-range antiferromagnetic order in the Hubbard model on a square 
lattice in the thermodynamic limit A 2 —> oo at hole concentration x = 0, temperature T = 0, and U/At > 0. Indeed, 
there is a large consensus that in the thermodynamic limit A 2 — > oo long-range antiferromagnetic order occurs in the 
ground state of the related isotropic spin-1/2 Heisenberg model on the square lattice^ 3 . - — . This implies that in the 
thermodynamic limit A 2 — > oo a similar long-range order sets in in the ground state of the half-filled Hubbard model 
on the square lattice at large U/At ^> 1 values. Moreover, for U/At > a similar ground-state order occurs in that 
limit in the ground state of the latter model, in terms of the spins of the rotated electrons. That according to the 
above property I the x — and m — absolute ground state is in the limit A 2 — ► oo invariant under the electron - 
rotated-electron unitary transformation then implies that in the thermodynamic limit TV 2 — > oo and for U/At > a 
long-range antiferromagnetic order sets in in the ground state of the half-filled Hubbard model on the square lattice 
in terms of the spins of the electrons as well. This agrees with many previous studies of that model, as for instance 
those of Refsi 11 ' 26 " — . As discussed below, there is strong evidence that both for x — 0, T > 0, and U/At > and 
for < x <C 1, T > 0, and U/At > such a ground-state long-range order is replaced by a ground-state short-range 
spiral-incommensurate spin order. 

B. x — and m = ground-state symmetry for N% —>■ oo and a necessary condition for its spontaneous 

symmetry breaking 

Our above arguments provide strong evidence that provided that in the thermodynamic limit iV 2 — > oo a long-range 
antiferromagnetic order occurs in the ground state of the related isotropic spin-1/2 Heisenberg model on the square 
lattice, a similar long range order sets in in that limit in the ground state of the half-filled Hubbard model on the 
square lattice for U /At > 0. However, our goal is not providing a mathematical proof that in the thermodynamic limit 
N% — > oo long-range antiferromagnetic order occurs in the half-filled Hubbard model on the square lattice. Although 
there is no such a proof, there is a large consensus that it should be so 11 ' 26 " — . An appropriate measure of such a 
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long-range order is the square of the staggered magnetization, 

(*GS\M S \*GS) = (*GS\ W £ W l*GS> = ^ £ C i ■ ( 3 °) 

\ a 3=0 J a j=0 

Here the spin operator sp. refers to the spin of an electron at the site of real-space coordinate fj in the original lattice, 
Cj is the spin correlation function in the spin-singlet S s = and x — ground state \%>gs) of the model for N% 3> 1 
very large but finite, 

Cj = (*Gs|*ft-»iil*Gs). ( 31 ) 

and = +1 if j refers to the sub-lattice A. Otherwise ej = — 1. This just compensates the sign of the antiferromagnetic 
correlation function (|3"TT) . For N% ^> 1 finite and even the x = and m = ground state |^gs) appearing in Eqs. 
(1501) and (|3"Tj) has zero momentum. 

The spin operator sp. appearing in Eqs. ([50]) and (|3"Tj) has operator Cartesian components s^ 1 , s^ 2 , and s~ 3 and 
refers to the spin of an electron at the site of real-space coordinate fj . In turn, we denote the corresponding rotated 
spin or spinon operator by sp. and its operator Cartesian components by s% 1 , s% 2 , and s% 3 . It refers to the spin of a 

3 3 3 3 

rotated electron at the site of real-space coordinate fj whose components appear in Eqs. (J7| and ([9|. Our spinous 
are such rotated spins. 

Under the emergence of long-range antiferromagnetic order in the limit N% — > oo, the square of the staggered 

magnetization (^gs\M s \^gs) extrapolates to a finite asymptotic absolute value \Coo\ of the correlation function 
(EH), 

lim (* G5 |^.|*Gfir) = |C 00 |. (32) 

ATJ->oo 

The related magnetic structure factor S(k) is the Fourier transform of that correlation function, 

JV a 2 -l 

S(k)= £ e^'Cj. (33) 

It can be measured directly in neutron scattering experiment o 3 ' 19 . As a result of its form (|33|) . S(if) = N%(M S ) will 
grow linearly with the number of sites N% if there is long-range antiferromagnetic order for N% — > oo, 

= (*gs\M s \*gs) = ^ + 0(l/N a ) . (34) 

Here the sub-lattice magnetization ttlaf = lmijv 2 ->oo 3(M S ) = 3|Coo| plays the role of antiferromagnetic order param- 
eter. 

The large scale DQMC calculations of Rebii provide useful information on the effective bandwidth, momentum 
distribution, and magnetic correlations of the half-filled Hubbard model on the square lattice. They employ the 
DQMC method, which provides an approximation-free solution of the such a model on square lattices large enough 
to use finite-size scaling to, for example, reliably extract the sub-lattice magnetization itiaf = 3|Coo| as a function 
of U /At. Such Monte Carlo calculations as well the random-phase approximation results of Refill reveal that ttlaf 
vanishes for If/At — > and is an increasing function of U/At that for approximately U/At ~ 2 saturates to the value 
m^f? ~ 0.614 of the isotropic spin-1/2 Heisenberg model on the square lattice determined in Refj2£. The quantum 
Monte Carlo DQMC results of Ref.— are an improvement of the corresponding results of Monte Carlo simulations 
of Ref. 26 , which predicted a lower value for the saturated itiaf- Below we relate the sub- lattice magnetization tuaf 
to an energy scale that can as well be used as antiferromagnetic order parameter. Combining that relation with the 
behavior tjiaf — in^p ps 0.614 for U /At 3> 1, we find the following approximate limiting behaviors, 



m A F ~ jj- e "v u ? U/At < 1 
hm 1(2-U/At\ 2 



lAF 4 

MM 



(^) ■ 



2f , U/At>2. (35) 
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where 

+ - 0.6142. (36) 

The expression given in Eq. (|55|) for u < U /At < 2 is also valid for < (u — U /At) -C 1 where u is a C/ /4t value 
found below to read uq « 1.302. Alike in Ref.— , Eqs. (|35|) and (|36|) refer to units where the classical Neel state has 
ttiaf = 1- In Ref^ units are used where that state has ttiaf = 1/2. In the latter units the parameter (|3f)|) reads 
instead ps 0.3071, consistently with the results of that reference. 

When expressed in terms of rotated-electron creation and annihilation operators, the Hamiltonian of the Hubbard 
model on the square lattice ([IJ has an infinite number of terms, as given in Eq. ([5]). For the half- filled model in the 
one- and two-electron subspace with both no rotated-electron double occupancy and no rotated-hole double occupancy 
the c fermion band is full for the ground state. Moreover, the excitations involving the emergence of c fermion holes 
are gapped. In the one- and two-electron subspace the model Hamiltonian ([5]) may then be expressed only in terms 
of spinon operators sp j whose operator components s l p , are given in Eq. ([?])• Consistently with the related results of 

Refi 7 -, one finds that up to fifth order in t/U the Hamiltonian ([5]) has in terms of such spin operators the following 
form, 

t 2 2t 4 
H = — > (sft. ■ Sff. — 1) — > fsi?. • sp. — 1) 

2JJ ' ' v r n r J2 > JJ3 / J v r 31 T J2 1 

(hh) (jih) 

t 4 

+ 21P 51 DjxjzDfoj^af^ ■ s? j3 - 1) 

31 ,32 ,33 

t 4 

— Sf H ■ Sj> ja — sp 3i ■ sp u — sV J2 • sp ja — Sf j2 ■ Sf u — sV 33 • Sf H ) 
5t 4 

jl ,32 ,33,3a 

+ (flf^j ' S rj i ){ s P j2 ' S f f 3 ' 3 ) — ( S f f 31 ' ^H 3 )(Sr j2 ' *Vj 4 )] • (37) 

Here (J1J2) refers to a summation running over nearest-neighboring sites and -Dj,,-/ = 1 for the real-space coordinates 
fj and fji corresponding to nearest-neigboring sites and Dj ji = otherwise. Analysis of the interactions in spin 
space of the Hamiltonian (|37|) reveals that some of its terms do not introduce frustration whereas other do. However, 
at half filling the spin interactions of the Hamiltonian (f3"T| including those of all higher order contributions do not 
destroy the sub- lattice magnetization tjiaf- They merely destabilize the classical Neel state, lessening the sub-lattice 
magnetization from its classical magnitude ttiaf = 1- 

That as obtained by different authors and method s 11 ' 23 ' 26 ' 27 the sub-lattice magnetization ttiaf of Eqs. ([34]) and 
(|35|) is indeed finite for the Hubbard model on the square lattice at U /At > provides strong evidence that for 
Nl -> 00 the spin correlation function (|31[) has long-range antiferromagnetic order. In contrast to ID, the quantum 
fluctuations associated with the interactions in spin space of the Hamiltonian (l37l) and its higher-order terms are 
not strong enough to destroy it. In turn, the exact Mermin- Wagner Theorem 3 ^ implies that at finite temperatures 
thermal fluctuations destroy such an order of the square-lattice model. 

A stronger confirmation is obtained by the scaling of the spectrum itself. It is an illustration of the mechanism 
of spontaneous symmetry breaking. Anderson was the first to point out that the spontaneous symmetry breaking 
mechanism that occurs in the thermodynamic limit ./V 2 — > 00 involves a whole tower of low-lying energy eigenstates of 
the finite system 37 . They collapse in that limit onto the ground state. One may investigate which energy eigenstates 
couple to the exact finite iV 2 >• 1 and x = and m = ground state |^gs) y i a the operator of the staggered 
magnetization M l s . Here I = ±, £3. We insert a complete set of energy eigenstates as follows, 

(* G 5|(A^i) 2 |^ GS ) = ^(^G^lA^il^X^lAfilxI/Gs) = l(*Gs|A^i|^)| 2 ; l = ±,x 3 . (38) 

i i 

Only excited energy eigenstates with momentum k — if and quantum numbers S v = 0, S c = N^/2, S s = 1, 
and S^ 3 = 0,±1 corresponding to / = 2:3, ± contribute to the sum of Eq. (|3"5)) . We emphasize that the quantum 
numbers S v — and S c = N%/2 remain unchanged and thus are the same as for the ground state I^gs)- We denote 
by \^>it) the lowest S s = 1, S v — 0, S c — N%/2, and k = [n,n] spin-triplet state whose excitation energy behaves 
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as 1/N% for finite N% 3> 1. For the range U/At > uq « 1.3 of interest for the studies of RefJ^ the contribution from 
such a lowest spin triplet state is by far the largest: For instance for approximately U/At > 2 the matrix-element 
square K^gsI-^II^it)! 2 exhausts the sum in Eq. (|3"5)) by more than 98.7%. (This is the value found by exact 
diagonalization for the related spin- 1/2 Heisenberg model on the square lattice in Reb^i. As similar result is expected 
for approximately U/At > uq m 1.3.) 

The special properties with respect to the lattice symmetry group of the lowest energy eigenstates contributing to 
the linear Goldstone modes of the corresponding S s = 1 spin-wave spectrum reveal the space-symmetry breaking of 
the N% — > oo ground state. In the present case of the half-filled Hubbard model on the square lattice the translation 
symmetry is broken. Hence both the k — [0,0] and k — [tt, tt] momenta appear among the lowest energy eigenstates 
contributing to the linear Goldstone modes of the S s = 1 spin-wave spectrum. 

The c and si fermion description can be used to derive such a spin-wave spectrum. For x = and m = it is 
generated in RefJ^ in terms of simple two- sl-fermion- hole processes. The investigations of that reference confirm 
that at x = the spin- wave spectrum includes linear Goldstone modes at momenta k = [0,0] and k = [tt, tt]. 
That the transverse spin-spin correlation function contains gapless poles, as predicted by the Goldstone Theorem, is 
consistent with in the limit N% — > oo the ground state breaking the continuous spin SU(2) rotational invariance of 
the Hamiltonian. The corresponding spin- wave spectrum is plotted in Fig. 1 of Rcf. 19 for the half- filled Hubbard 
model on the square lattice in the thermodynamic limit N% — > oo. (As mentioned above, for large N% the lowest 
k = [tt, 7r] excitation energy vanishes as 1/N^.) The occurrence of linear Goldstone modes at momenta k — [0, 0] and 
k = [7r,7r] in the theoretical spin-wave spectrum derived in RefJ^ explicitly confirms that the results of the c and si 
fermion description used in the studies of this paper are fully consistent with for x — 0, m = 0, and temperature 
T = long-range antiferromagnetic setting in as A^ y oo. 

Importantly, for U/At w 1.525 and t m 295 meV the spin-wave spectrum of the parent compound La2Cu04 (LCO)£ 
is quantitatively described by the corresponding theoretical spectrum derived in Ref. — in terms of simple spinon pair 
breaking si fermion processes. Within the present status of the scheme used in the studies of that reference one 
cannot calculate explicitly matrix elements of the two-electron spin-triplet operator between energy eigenstates and 
corresponding spectral- weight distributions. The x = and m — results of Ref^- on the spin- wave spectrum of the 
parent compound LCO profit from combination of the c and si fermion description with the complementary method 
of Rcf.^. They reveal that the microscopic mechanisms that generate the coherent spectral- weight spin- wave energy 
spectrum are in terms of spinon pair breaking si fermion processes very simple. Indeed the two-spinon si fermion 
description renders a complex many-electron problem involving summation of an infinite set of ladder diagrams^ into 
a non-interacting two-sl-fermion-hole spectrum, described by simple analytical expressions. 

Within the semi-classical description of spin waves, they can be pictured as long wave-length twists of the order 
parameter. In turn, within the present quantum description a spinon bond pair of a spin-singlet two-spinon si fermion 
is broken, giving rise to two independent spin-up spinons or two independent spin-down spinons. All remaining 
spinons in the problem remain confined with spin-neutral two-spinon si fermions. The two deconfined spinons are 
invariant under the electron - rotated-electron unitary transformation. Thus their spin-triplet excite-state occupancies 
correspond to an isolated vanishing- energy and vanishing-momentum mode below a continuum of two-sl-fermion-holc 
excitations. Indeed, the momenta ±<f corresponding to the spinon relative motion in the spinon pair of the broken si 
fermion are transferred over to two holes, respectively, that emerge in the si band. Under spin-triplet excitations such 
a two-spinon si fermion breaking may introduce an extra momentum contribution Sq. It corresponds to the motion 
of the center of mass of the spinon broken pair. In the latter case the two emerging si band holes have momenta q' 
and q" given in Eqs. (gH) and (|2T|) . 

As found in Ref. 1 — , the processes associated with most momenta q' and q" lead to an incoherent background of 
spin spectral weight. In turn, the spin coherent spectral- weight distribution refers to the spin- wave spectrum. The 
processes that generate such a coherent spin weight are such that one of the si band momenta q' and q" belongs to 
the boundary line and the other points in a nodal direction. That some of the spin spectral weight is incoherent is 
consistent with the sub- lattice magnetization obeying the inequality wiaf < 1j rather than reading rriAF = 1, as for 
the classical Neel state. In turn, we argue below that the spin spectral weight is fully incoherent for spin excitations 
of m = and x > ground states. 

In the thermodynamic limit N£ — > oo a large number of low-lying energy eigenstates \^i) with momentum k = ir 
and quantum numbers S v = 0, S c = N%/2, S s = 1, and Sf 3 = 0, ±1 that contribute to the sum of Eq. (|3"5)) converge 
to the corresponding N% 3> 1 finite ground state. To illustrate the mechanism of spontaneous symmetry breaking, we 
consider for simplicity that the lowest-energy spin triplet state I^it) belonging to that set of energy eigenstates gives 
rise to such a symmetry breaking. Indeed, for intermediate and large U /At values such a state exhausts the sum in Eq. 
(f3"5f by more than 98.7%. However, the lower broken symmetry of the final ground state is the same independently of 
the number of low-lying states considered in the analysis of the problem. In the presence of a small staggered field B s 
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in the ^-direction, a new ground state emerges due to the additional term M^ 3 B X:j in the Hamiltonian. This state, 

|**-gs) « O g <?|*gs) + Cit|*it) , (39) 

has a finite staggered magnetization, 

M s = (y b _ GS \M:\y b ~Gs) = 2C G sC 1t ^>gs\M z s \^it) ■ (40) 

In turn, (* GS | M S Z |* GS ) = and {V vr\M%\$ vr) = 0. The new ground state |* 6 _ G s) has not a well defined spin 5 S , 
yet it has the same quantum numbers S n = and S c — N%/2 as the spin-singlet ground state |*gs)- 

In the thermodynamic limit N£ —> oo, the system assumes the largest possible magnetization for arbitrarily small 
staggering field B s , so that C GS = C XT = I/a/2 and thus M s = (* GS |M s z |*i T ). The state \^ b -Gs) of Eq. ^ 
is contained in a reduced subspace spanned by |$gs) an d |$ir). In such a reduced subspace the unity partition 
operator reads s» |* G s)(vE' G s| + I^itX^itI- We then obtain the following staggered magnetization squared, 

M 2 = ^ b ^ GS \M z s \^ b -GS? = (^gs\(M!) 2 \^gs) ■ (41) 

It follows that the staggered magnetization squared in the iV 2 — > oo ground state |\&&_gs) with broken symmetry is 
identical to the long-range antiferromagnetic order of the correlation function in the spin-singlet N% 3> 1 ground state 
I^Gs), consistently with Eq. (1521 . Indeed, the studies of Ref. 12 reveal that for a number of sites iV 2 ^> 1 large but 
finite the m = ground states are spin-singlet states. For the present case of a m — and x — ground state this 
agrees with an exact theorem introduced and proved in Ref£. 

The new ground state \^b-Gs) of Eq. (fBTJf is a superposition of states with different spin S s and thus breaks 
the Hamiltonian spin SU(2) symmetry contained in its global 5*0(3) x 50(3) x [7(1) = [SU(2) x SU(2) x U(l)]/Z% 
symmetry. Specifically, the spin rotational symmetry SU(2) is spontaneously broken to [7(1) by the formation of the 
staggered magnetization. In turn, the new ground state has the same quantum numbers S v = and 5 C = N%/2 as the 
spin-singlet ground state I^gs)- Therefore, the corresponding ?7-spin symmetry SU(2) and c fermion symmetry [7(1), 
respectively, are not broken. A similar result is obtained if besides \^gs) an d l^ir), the new ground state \^b-Gs) 
contains a larger set of low-lying energy eigenstates with momentum k = 7? and quantum numbers 5 r; = 0, 
5 C = Nl/2, S s = 1, and 5f = 0,±1 other than |* 1T ). Hence rather than SO (3) x 50(3) x 17(1), the symmetry of 
new ground state |* 6 -gs) is [[7(2) x U(l)]/Z$ = [50(3) x [7(1) x U{l)]/Z 2 = [SU{2) x [7(1) x [7(1)]/Z 2 . 

Consistently with the studies of the x = two-spinon si fermion pairing energy of Ref.—, in the Hubbard model 
on the square lattice state configurations such that iV 2 ^ = [TV 2 — 25 c ] = and thus N% a = -/V 2 exist below an energy 
scale n° = lim^^o Such an energy scale refers to the spin degrees of freedom yet it equals one half the charge Mott- 
Hubbard gap. The latter gap defines the range /i g (—A* , A* ) of the chemical potential /i at x = and m = 12 i 19 . The 
energy scale /jP may be used as order parameter of the long-range antiferromagnetic order. It equals the excitation 
energy below which such an order exists in the limit TV 2 — > oo at m = 0, x = 0, and zero temperature T = 0. The 
limiting behaviors of /jP are approximately the following, 

fjP 32 te^vf , [7/4* < 1, 

« — Vl + (£7/4*-u ), wo < ?7/4* < m , 

7T 

« [[7/2 - At] ; [7/4* > 1 . (42) 

It vanishes /jP — > for [7/4* — > and is finite and an increasing function of U /At for [7/4* finite. For [7 /4f — > oo it 
behaves as fi° s» [7/2 — > oo. This is why in that limit, when the spin degrees of freedom of the half-filled Hubbard 
model are described by the isotropic spin- 1/2 Heisenberg model, the state configurations for which N% s = -/V 2 exist 
at any finite energy. 

The behavior /j, (U/At) « A*o( M o)\/l + (U /At — uq) reported in Eq. (|42|) is expected to be a good approximation 
for the intermediate [7/4* range [7/4* G (uoi u i) of interest for the square-lattice quantum liquid studies of Refi2. 
Here uo ~ 1.302 and u\ ~ 1.600. The approximate magnitude Ato(wo) ~ [2e 1 /7r] * is that consistent with the relation 
A t o('"o) ~ A*o(' u *)/'\/l + (u* — u o)- The value [7/4* = = 1.525 is that at which the studies of that reference 
lead by a completely different method to ^o(w*) ~ 566 mcV for * = 295 meV and U/At s» u* = 1.525. The use of 
Mo(wo) ~ {2e 1 /ir}t in the formula ^ Q (U/At) w Mo("o) V 1 + (^/ 4i _ u o) leads for * = 295 meV, U/At ss = 1.525, 
and ito = 1.302 to nearly the same magnitude, /zo(m*) ~ 565 meV. In turn, the value [7/4* = uq — 1.302 is that at 
which the energy parameter 2A is found below to reach its maximum magnitude. That energy parameter such that 
2Aq < fj,Q for U/At > is the energy below which the short-range incommensurate-spiral spin order considered below 
survives for < x <C 1, m = 0, and zero temperature T = 0. 
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We make the reasonable assumption that at T = the energy parameter fjP and the sub-lattice magnetization tuaf 
of Eq. (|35|) scale in the same way as follows, 



fJ }' = Um A Fa". (43) 

A naive spin-density wave mean- field approach^ 7 , leads to the relation fi° — [U/2] rriAF- In turn, in the case of the 
relation (j43|) the coefficient a° = fi°/[U tuaf] has the following limiting behaviors, 

a a = 0.500 , U/At < 1 , 

pa 0.603, U/At = u a Pa 1.302, 

0.536, t//4i = u* ps 1.525, 

pa 0.519, U/At = u t w 1.600, 

= « 0.814, [//4t»l. (44) 

AF 

The physical quantities m^F of Eq. (j35"|) , Uttiaf, and /x° of Eq. (|42|) are increasing functions of U/At. In turn, 
a = /jP/Uttiaf is a monotonous function of U/At with both minima and maxima. In the limit U/At — » the effects 
of all fluctuations vanish and the value a — 1/2 is that predicted by mean-field theory relation ijP = [U/2] tjiaf- As 
a function of U/At, the coefficient a first increases until reaching a maximum at a f/ /4t value below uo- Interestingly, 
in the intermediate range U/At € (uq, u±) it is a decreasing function of U/At. It reaches a minimum value larger than 
0.500 and smaller than 0.519 at a U/At magnitude laying between U/At = u\ « 1.6 and U/At = 2. It then becomes 
an increasing function of U/At, reaching the value a pa 0.814 for U/At 3> 1. Since /jP/U —> 1/2 as U/At — > oo, note 
that the corresponding U/At ^> 1 expression a = 1/[2itiaf] would give a = 1/2 if the tuaf value was that of the 
classical Neel state tuaf = 1, rather than that of the spin- 1/2 Heisenberg model tuaf = m AF ~ 0.6142. 

Consistently with for the Hubbard model on a square lattice [i° being the energy below which the long-range 
antiferromagnetic order exists in the limit N% — »ooata; = 0,m = 0, and zero temperature T = 0, the limiting 
behaviors /x° pa 32 i e - 27r \/*A 7 an( j „ jj j 2 

given in Eq. (j42]) for ?7/4t -C 1 and U/At » 1, respectively, are those 
of the zero-temperature spin gap of Eq. (13) of Refi^. Within the operator description used in this paper the c 
fermions and spinon occupancies generate the state representations of the groups U{\) and spin SU(2), respectively, 
in the model global 50(3) x 50(3) x 17(1) = [SU(2) x SU(2) x C/(1)]/Zf symmetry. In turn, the studies of Ref> 
isolate strongly fluctuating modes generated by the Hamiltonian (JXJ) Hubbard term according to the charge i7(l) and 
spin SU{2) symmetries. Within the gauge transformation of that reference, the strongly correlated problem is casted 
into a system of noninteracting u h fermions" submerged in the bath of strongly fluctuating U(l) and SU(2) gauge 
potentials. Those couple to fermions via hopping term plus the Zeeman-type contribution with a massive field g(fj-r). 
Within the description of that reference a U(l) and SU(2) gauge transformation is used to factorize the charge and 
spin contribution to the original electron operator in terms of the emergent gauge fields. The U(l) charge h fermions 
and SU{2) spins of Ref.— refer to the c fermions and spinons, respectively, of our description. In what the relation 
of the energy scale of Eq. (J42J to the results of that reference is concerned, in the x = antiferromagnetic phase the 
above massive field gifjr) assumes the staggered form, 

g(f j T) = n e iS ^ . (45) 

The description of Refi 3 ^ is valid for large £7/4£ values, so that /io ~ U/2, as given in Eq. (|42|) . In turn, the c 
fermions and spinons of this paper are directly related to the rotated electrons of the description of Ref J^, whose 
double occupancy is a good quantum number for U/At > 0. This is why the energy scale /io of Eq. (|45|) is here well 

defined for the whole range of U/At values. Consistently, fjP pa 32 te~^^v —> and thus p(fj t) for U/At —> 
upon the disappearance of the long-range antiferromagnetic order, as given in Eq. (|42l) . 



C. Short-range incommensurate-spiral spin order at x = and < T <C Tq and for < x 1 and < T <C To 

Here it is argued that the spin effective lattice being identical to the original lattice is a necessary condition for the 
emergence of a ground-state long-range antiferromagnetic order in the limit N% — > oo. Moreover, it follows from the 
property II of Section III-A that since the ground-state rotated-electron occupancy configurations are for finite values 
of U/At more ordered than those of the electrons generating the same state, a lack of long-range antiferromagnetic 
order of the rotated-electron spins in the limit N% — > oo for N% /N% < 1 would imply a similar lack of long-range 
antiferromagnetic order for the spins of the original electrons as N% — > oo. The inequality N% /N% < 1 applies to the 
x > and m = ground states. For those the spin effective lattice is different from the original lattice, in contrast to 
that of the absolute x = and m = ground state. 
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1. The general effects of hole doping accounting for the hidden global (7(1) symmetry 

Quantum-Monte Carlo methods, when applicable, are the only unbiased tools for quantitative studies of the effects 
of a small hole concentration x in the physics of the Hubbard model on the square lattice. Unfortunately, some of 
the interactions that become active for x > cannot be studied by Quantum-Monte Carlo simulations due to the 
well-known "sign problem" . While part of our results are argued on phenomenological grounds, taking account for 
the effects of the hidden global U(l) symmetry found recently for the Hubbard model on any bipartite lattice in Rcf. 4 
introduces a new scenario and framework, which may be useful for future quantitative studies of the square-lattice 
model. In addition, such a new scenario allows the preliminary qualitative discussion of the problem presented in this 
paper. 

In accordance to a general theorem proved in Ref.— , at half filling the terms of the Hamiltonian ([5]) expansion in 
t/U with odd powers in t vanish due to the particle- hole symmetry and the resulting invariance of the spectrum under 
t — > —t. For instance, the x = and m — Hamiltonian terms of Eq. (157)) result only from the terms of order t 2 
and t A of the Hamiltonian (6) of Ref^. (The To and T±\ operators of that reference include a factor t absent in the 
corresponding operators of Eq. ©.) In turn, for finite hole concentration x > and vanishing spin density m = the 
expansion in powers of t/U of the Hamiltonian ([S]) involves terms with odd powers in t, absent at x = 0. We argue 
that the emergence of such new terms absent at x = along with related effects associated with changes in the spin 
effective lattice upon "turning on" the hole concentration x destroy the sub-lattice magnetization tuaf- 

Expression of the Hamiltonian of the Hubbard model on the square lattice (|TJ) in terms of rotated-electron creation 
and annihilation operators leads to the Hamiltonian (JS|), which has an infinite number of terms. Its expansion up 
to fifth order in t/U leads for the x — half-filled Hubbard model in the one- and two-electron subspace to the 
Hamiltonian expression given in Eq. (|37p . which involves only spinon operators. In turn, for x > the Hamiltonian 
([5]) may be expressed in terms of c fermion operators, 77-spinon operators, and spinon operators. This is done by 
the use of the expressions provided in Eq. (|A1[) of Appendix A for the rotated-electron operators in terms of such 
operators. The uniquely obtained Hamiltonian also has an infinite number of terms and thus is rather complex. 

The problem slightly simplifies for the x > Hamiltonian ([5]) in the one- and two-electron subspace. In terms 
of rotated-electron creation and annihilation operators, its terms generated up to fourth order in t/U are for x > 
and within a unitary transformation the equivalent to the t — J model with ring exchange and various correlated 
hoppings 13 . Moreover, in the one- and two-electron subspace such an Hamiltonian may be expressed in terms of only 
c fermion operators and spinon operators. The rotated quasi-spin operators are given by q l p. — s|» +p l p. where I = ±, z. 

Here s^. and p l ~. are the spinon local operators and 77-spinon local operators, respectively, given in Eq. ([7j. Since in 
the one- and two-electron subspace there are no rotated-electron doubly occupied sites, the 77-spin effective lattice is 
empty and the operator p l ~_ plays no active role. However, the obtained Hamiltonian expression is for x > much 
more involved than that given in Eq. (|37[) for x = and is omitted here. In addition to involving spinon operators, it 
contains c fermion operators and its number of terms is much larger than for x = 0. Moreover, it involves c fermion 
- spinon interactions and both non-frustated isotropic spinon interactions and spinon interactions with some degree 
of frustration. Hence such interactions cannot be studied by Quantum-Monte Carlo simulations due to the "sign 
problem" . 

In such an involved x > Hamiltonian expression the c fermion operators act onto the c effective lattice occupancies 
and the spinon operators act onto the spin effective lattice occupancies. (The c fermion operators commute with the 
spinon operators.) Also in the x > Hamiltonian (JSJ) general expression containing an infinite number of terms, the 
c fermion operators act onto the c effective lattice occupancies, the 77-spinon operators act onto the 77-spin effective 
lattice occupancies, and the spinon operators act onto the spin effective lattice occupancies. The state representations 
of the new hidden global U(l) symmetry are generated by the c fermion occupancy configurations in the c effective 
lattice. As mentioned above, the point is that for each energy eigenstate the relative positions of the sites of the 77-spin 
and spin effective lattices in the original lattice are stored in the c effective lattice: For U /At > the latter lattice 
is identical to the original lattice for all A N °- energy eigenstates and the N c = 2S C c fermions occupy the 2S C sites 
singly occupied in the original lattice by the rotated electrons of the description introduced in RefJ^. We recall that 
for U /At > the sites of (i) the 77-spin effective lattice and (ii) the spin effective lattice correspond to those (i) singly 
occupied and (ii) unoccupied and doubly occupied by rotated electrons. Hence for each energy eigenstate the relative 
positions of the sites of the (i) 77-spin effective lattice and (ii) spin effective lattice in the original lattice refer to the 
relative positions of the sites of the c effective lattice (i) unoccupied by c fermions and (ii) occupied by c fermions, 
respectively. 

That the relative positions of the sites of the 77-spin and spin effective lattices in the original lattice are stored in the c 
fermion occupancy configurations of the c effective lattice is consistent with for N% — > 00, N^/N 2 = [1 — 2S C /N%] > 0, 
and N^JNa = 2S C /N 2 > the occupancy configurations of the 77-spinons and spinons of the general description of 
Rcf. 12 referring to independent //-spin and spin effective lattices, respectively. Moreover, in that limit these lattices 
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can be considered as square lattices with spacing a v and a s , respectively, provided in Eq. (|A8|) of Appendix A. For the 
one- and two-electron subspace the latter spacing reads a s — a/y/1 — x, as given in Eq. (jl8p . Hence alike for x = 
and m — ground state, for x > and m — ground states the the M s = 2S C spin-1/2 spinons remain occupying a 
full spin effective lattice with N% = M s = 2S C sites. 

Although in the limit N% — > oo the occupancies of the spin, ?7-spin, and c effective lattices are independent, for 
the x > Hubbard model on the square lattice in the one- and two-electron subspace the M s = 2S C spin-1/2 
spinons interact with each other and with the c fermions. However, as mentioned above the expression of the x > 
Hamiltonian ([5]) in the one- and two-electron subspace in terms of c fermion operators and spinon operators leads to 
a very involved quantum problem. Indeed, the usefulness of the c, spin, and si effective lattices and corresponding c 
and si momentum bands description refers to that Hamiltonian in the one- and two-electron subspace in normal order 
relative to the initial x > and m = ground state. Fortunately, such a ground-state normal-ordering simplifies the 
quantum problem^. It provides implicitly and naturally a criterion for the selection of a few dominant Hamiltonian 
terms expressed in terms of c fermion operators and two-spinon si fermion operators. That problem is studied in 
RefJ^ in terms of a suitable energy functional valid for intermediate and large values of U/At. The spin degrees of 
freedom of such a functional describe both the x = and x > problems. The model spin spectrum relative to the 
x = and m — ground state is one of the few problems for which there are results from controlled approximations 
involving summation of an infinite set of ladder diagrams^. As mentioned above, the spin spectrum provided by 
the energy functional of RefJ^ quantitatively agrees with both the spin-wave spectrum derived in Refi^£ and that 
observed in the parent compound LCO 3 . 

Alike for the x — and m — ground state, for x > and m — ground states the M s = 2S C spinons 
remain confined within N s \ = S c spin-neutral two-spinon si fermions. The si fermions are generated from bosonic 
spin- neutral two-spinon si bond particles and thus have long-range interactions associated with the effective vector 
potential A s i(fj) of Eq. (|A2I) of Appendix A. An important property is that for x > and m = ground states 
the si fermion momentum band is full and for one-electron and two-electron excited states displays a single hole and 
none or two holes, respectively-^. The si - si fermion interactions associated with the effective vector potential of Eq. 
(IA2[) of Appendix A are stronger than those that arise between the emerging si fermions and pre-existing c fermions. 
In spite of that, the former do not lead to si - si fermion inelastic scattering. The obvious reason is that due to 
phase-space restrictions associated with the exclusion principle and energy and momentum conservation requirements 
there are no available momentum values in the si band for excited-state occupancy configurations. 

Both for the x = and m = ground state and x > and m — ground states the M s = 2S C spin-1/2 spinons 
occupy a full spin effective lattice with N% = M s = 2S C sites. However, an important point is that within the present 
description there is a qualitative difference between the x = and x > problems. At x = the spin effective lattice 
is identical to the original lattice whereas for x > it has a smaller number of sites N% = 2S C < N% so that its 
spacing a s = a/y/1 — x is larger. Within the description of the problem used in the studies of this paper, this is 
one of the main effects of hole doping. Indeed, the lack of the 77-spin effective lattice, N% = x N% = 0, occurring at 
x = implies that the spin effective lattice is identical to the original lattice and thus has the same number of sites 
N% = 2S C — N% as that lattice. We argue that the spin effective lattice being identical to the original lattice is a 
necessary condition for the ground-state long-range antiferromagnetic order to emerge as N% — > 00. The concepts 
of spin effective lattice and si effective lattice are only valid in that limit. Only in it are such lattices approximate 
square lattices with spacing a s = a/y/1 — x given in Eq. (|18[) and a s \ ~ y/2a s provided in Eq. (f2~Tj) . respectively. The 
incommensurability relative to the original square lattice spacing a of the effective lattice spacings a s = a/y/1 — x 
and a s i ~ V2a s is consistent with a ground-state long-range antiferromagnetic order occurring for N% — > 00 only at 
x = 0. Such an incommensurability allows processes that destroy long-range antiferromagnetic order becoming active 
for x > 0. Within the expansion in powers of t/U of the Hamiltonian ([5]), such processes are associated for instance 
with Hamiltonian terms with odd powers in t, absent at x = 0. 

The possible existence of low-lying spin excitations for x > is not a sufficient condition for a spontaneously broken 
symmetry to occur in the limit N% — > 00. According to the results of RefJ^, the m = and x > ground states 
are for N% » 1 large but finite spin-singlet states. Note though that only if the state obtained by application of the 
two-electron spin-flip operator onto such ground states had finite overlap with low-lying spin excitations would these 
states acquire a long-range spin order as N% — > 00. However, we argue that if low- lying spin excitations exist for 
m = and x > such an overlap vanishes. It vanishes as well in the case of non existence of low-lying spin excitations 
due to a spin gap: There is no coherent spin spectral weight both for vanishing and/or finite energy. Therefore, 
in addition to the m — and x > ground states remaining spin-singlet states in that limit, the corresponding 
spin-triplet spectrum is fully incoherent: Its sharp spectral features are not <5-function like. The absence of coherent 
spin- wave excitations is consistent with the lack of long-range spin order in the m = and x > ground states in the 
limit N% -t 00. 

For a x > and m = ground state one has that the number of sites of the 77-spin effective lattice N% = x N% > 
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is non zero even if the hole concentration < x -C 1 is very small. For such a state the spin effective lattice has 
a number of sites N% < iV 2 smaller than that of the original lattice. Consistently with the corresponding effective 
lattice spacings reading a s = a/y/1 — x and a s i ~ V2a s , below further evidence is provided that for < x -C 1 the 
ground state of the Hubbard model on the square lattice has a short-range incomensurate-spiral spin order. Moreover, 
according to the results of Ref.*— the ground state has a short-range spin order for < x < x*. For x > x* it is a 
disordered state without short-range spin order. Here x* > 0.23 for approximately U/4t > u ~ 1.3. 

Finally, we provide further evidence that the form of the si effective lattice spacing (|2~Tj) is for the Hubbard model 
on the square-lattice consistent with the above mentioned x = and x > spin orders, respectively. That at x = 
and m — the spacing of the square si effective lattice is given by a s \ = \/2a reveals that then its periodicity has 
increased relative to that of the original lattice, which in that case is identical to the spin effective lattice. Indeed, 
at x — the si effective lattice is one of the two sub-lattices of the original lattice and thus refers to a \/2 x \J~2 
reconstruction in which the periodicity of the spin-sub-system real-space structure is increased. Such an effect is 
consistent with the occurrence of the long-range antiferromagnetic order for N% — > oo at x = and m — 0. In turn, 
that for x > and m = the square si effective lattice spacing reads instead a s i ~ •v/2/(l — x) a is consistent with 
the emergence of the short-range incommensurate-spiral spin order. Indeed now the si effective lattice is one of the 
two sub-lattices of the spin effective lattice which for x > is different from and incommensurate to the original 
lattice. 



2. Quantum, and thermal phase transitions 

It is argued above that the lack of finite overlap of the state generated by application of the two-electron spin- 
flip operator onto a x > and m = ground state with low-lying excited states is behind such a ground state 
remaining a spin-singlet state for N% — > oo, alike for TV 2 ^> 1 large but finite^. This is so independently of the 
existence or non existence (spin gap) of spin-triplet low-lying states. Hence such a ground state has no long-range 
antiferromagnetic order in the limit N% — > oo and its symmetry is that of the Hamiltonian, SO(3) x 50(3) x £7(1) = 
[SU(2) x SU(2) x U(X)\jZ\. A necessary condition for the occurrence of such a long-range order is according to our 
above analysis that the spin effective lattice is identical to the original lattice. That condition is not met by x > 
and m = ground states. For them the Hamiltonian hidden global U(l) symmetry generator eigenvalue S c obeys the 
inequality S c < N%/2. It then follows that the number of sites of the 77-spin effective lattice N% = [N% — 2S C ] > 
is finite, so that the number of sites of the spin effective lattice N% = 2S C < N% is smaller than the number of 
sites of the original lattice. For the one- and two-electron subspace considered in this paper such inequalities read 
N% = x Nl > and /V 2 . = (1 - x) N% < N%, respectively. 

There is strong evidence of the occurrence in the half-filled Hubbard model on the square lattice of strong short- 
range antiferromagnetic correlations for finite temperatures T > below a crossover temperature called T x in Refi^, 
which here we denote by Tq . This is consistent with then the system being driven into a phase with short-range 
spin order. Furthermore, that the occurrence of long-range antiferromagnetic order as N% — > 00 requires that T = 0, 
N% = 0, and N% = N% is consistent with the short-range spin order occurring for m = 0, < x <C 1, and < T < T * 
having basic similarities to that occurring for m = 0, x = 0, and < T < T *. The latter order was studied previously 
in Ref. 29 for < T < T *. 

As further justified below, for both vanishing and finite temperatures a phase displaying a short-range spiral- 
incommensurate spin order is then expected to occur for (i) m = 0, < x <C 1, and < T « T ' and (ii) m = 0, 
x = 0, and O-cT^cTq. Atm = and temperatures below T *, the system is driven both for (i) < x <§; 1 and 
< T < Tq and (ii) x = and < T < T * into a renormalized classical regime where the N% — > 00, x — 0, and 
T = long-range antiferromagnetic order is replaced by such a short-range spin order, which is a quasi-long-range 
spin order as that studied in Ref4^ for simpler spin systems. 

An interesting physical issue is whether the quantum phase transition separating the x = and m = ground state 
from the < x <^ 1 ground state at T = corresponds to a "deconfined" quantum critical point3. Indeed, both at 
x = and for < x <Si 1 the ground-state M s — 2S C spinons are confined within A^ s i — S c spin-neutral two-spinon 
si fermions. When expressed in terms of rotated-electron creation and annihilation operators, the Hamiltonian of 
the Hubbard model on the square lattice has for U /At finite an infinite number of terms, as given in Eq. ([5]). As 
mentioned above, the Hamiltonian terms generated up to fourth order in t/U are for x > and within a unitary 
transformation the equivalent to the t — J model with ring exchange and various correlated hoppings*-. In addition 
to non-frustated isotropic spin interactions as those considered for a spin- 1/2 Heisenberg model on the square lattice 
in Refsi 21 ! 22 , the present spinons have some degree of frustration as well. Hence their interactions cannot be studied 
by Quantum-Monte Carlo simulations due to the "sign problem" . A related interesting open question is whether for 
x > and m = the the spin degrees of freedom of the present square-lattice quantum liquid refer to a valence-bond 
solid or some type of related valence-bond liquid. In either case, that both at x = and for x > the M s = 2S C 
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spinons are confined within N s i = S c spin-neutral two-spinon si fermions strongly suggests that the corresponding 
quantum phase transition refers indeed to a "deconfined" quantum critical point. If this is so, the critical point is 
characterized by deconfined spin-1/2 spinons coupled to some emergent £7(1) gauge fielder— . 

We denote by 2Ao the energy below which the short-range incommensurate-spiral spin order with strong antifcr- 
romagnetic correlations survives at zero-temperature, m = 0, and < i C 1. That energy parameter has a U/At 
dependence qualitatively similar to that of the energy scale 2/c^Tq*, with the equality 2Ao ~ 2fcsTg approximately 
holding. Except for U/At — > 0, such an energy scale has a different origin than the order parameter fjP = U tuaf oP 
of Eq. (|43l) . which is proportional to the sub- lattice magnetization uiaf of Eq. (|3"5"1) . Indeed, for < x <C 1 the lack 
of a long-range antiferromagnetic order implies that uiaf = 0. 



3. The U /At dependence of the 0<i<l energy scales 

The energy scale 2A ~ 2/cbTq considered here plays a major role in the square-lattice quantum-liquid studies 
of Ref. 19 . Here we address the problem of its £7/44 dependence by combining the results obtained from the use of 
our general description with those of the low-temperature approach to the half-filled Hubbard model on the square 
lattice of Refj2£. The investigations of that reference focus on temperatures < T -C Tq • The energy parameter 
2A « 2fcs T * refers to the limit 2A = lim^^o 2| A| of an x dependent energy scale 2| A| that plays the role of order 
parameter of the phase with short-range spin order. For < x <C 1 and intermediate and large U / At values such an 
energy scale reads, 

2|A| « 2A ^1 - , 0<x<l, £7/44 > u Q « 1.302 . (46) 

The linear dependence on x of [2A - 2|A|] w (x/xl) 2A for < x < 1 is justified in RefJ£. Here x® w 2r s /ir, the 
ratio r s = 2A /8M / ° 1 plays an important role in the square-lattice quantum liquid, and = lim^o W s i = W s i|aj=o 
where W s x is the nodal energy bandwidth W s \ of the si fermion dispersion defined in Ref.—. Its maximum magnitude 
is reached at U /At — 0. For £7 /At > it decreases monotonously for increasing values of £7/44, vanishing for £7 /At — > oo. 
That for U/At — > oo both W s \ — > and |A| — > is associated with the full degeneracy of the spin configurations 
reached in that limit. In it the spectrum of the two-spinon composite si fermions becomes dispersionless. The limiting 
behaviors of the m = energy parameter 8W^ = lim x _>.o 8W s i = 8W s i|a; = o contributing to the ratio r s = 2A /8W S ° 1 
readi 9 -, 

8W° X = 32t , U/At = ; %W° sl w | M. , U/8 2 t > 1 . (47) 

In contrast to the x > energy parameter 2|A|, the energy scale is well defined both at x = and for x > 0, 
having the same magnitude at x = and for x — > 0. Its magnitude = [49.6/295] t w 0.168t obtained 

at U/At — u* — 1.525 in Ref,— is about 12.25 times smaller than that found by use of the limiting expression 
Airt 2 /U w 2.060 1 of Eq. at U/At = 1.525. This reveals that such a limiting expression is valid for a smaller 

range of very large U/At values than in ID. Indeed, Ant 2 /U w 0.168 1 at U/At w 18.70, whereas Wjh w 0.168 1 at 
?7/4t w 1.525. Therefore, the relation = J w Ant 2 /U is valid for approximately t//8 2 < > 1, as given in Eq. ([17]). 
For the Hubbard model on the square lattice the energy scale J w Ant 2 /U controls the physics for a smaller [7/44 
range than in ID, which corresponds to very large U/At values such that U/8 2 t 3> 1. Hence the intermediate- U/At 
range plays a major role in the physics of that model, as confirmed by the related studies of RefJ^. 

That for U/At > the energy scale 2A = lim :E ^o2|A| has a different origin than the order parameter fi° = 
U rriAFoP of Eq. (|43[) is consistent with for U/At > their magnitudes being different. However, symmetry arguments 
related to the disappearance of the N 2 — > oo and x = ground-state long-range antiferromagnetic order for U/4t —> 
imply that limu/it-to 2Ao = (J* • Indeed the ratio 2Ao//i° — > 1 involving the two energy scales becomes one in 

that limit. Moreover, the energy scale 2A interpolates between 2A = /i° ~ 32t e^^V 4 */ 17 for U/At <C 1 and 
2Ao = SW^ w 7r [8t] 2 /U for approximately U /8 2 t 3> 1. It goes through a maximum magnitude at a [7/44 value found 
below to be approximately given by U/At — uq w 1.302. 

The energy parameter is an increasing function of U/At. As given in Eq. (|42l) . it behaves as //o ~ 32ie _7r v / ^" 
for {7/44 -C 1 and as /x ~ [£7 — 84] for [7/44 3> 1. In turn, the energy scale 8W^ is a decreasing function of L7/44. 
According to Eq. (gTJ it is given by SW^ = 324 at U/At = and decreases approximately as 8W° t 7r [84] 2 /L7 for 
very large U/At. Since the energy scale 2Ao interpolates between these two behaviors, it vanishes both for U/At ^ 

and [7/44 — > oo as 2Ao ~ 324 e~ 71 V 4 '' /c/ and 2A m n(8t) 2 /U, respectively. In these two limits it becomes fiQ and 
the energy scale 8W^ X ~ n(8t) 2 /U associated with the strong 0<i<1 antiferromagnetic correlations, respectively. 
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This is consistent with its maximum magnitude being reached at an intermediate U/At value U/At = uo at which 
the equality fio ra 814^ holds. We then define such a U/At value as that at which the lines 814^ = 8W^ 1 {U/At) and 
= V>o(U/4t) cross and thus the equality /j,o(uq) — 8W s ° 1 ('Uo) holds. The ratio r s — 2A /8Wg 1 is parametrized in 
the following as r s — e~ As where X s = | ln(2Ao/8M / s ° 1 )| controls such an interpolation behavior. 
The energy scale 2Aq can then be expressed as, 

2A = r s 8Wj x = 8W° ± e" A * ; A s = | ln(2A /8^ s ° 1 )| , (48) 

where A s has the limiting behaviors, 



A s = wy/At/U, U/U<1; \ s xi4tu a /U, u 00 <U/At < in; X s =0, U/At^oo, 
uoo w (wo/tt) 2 « 0.171 ; u ra 1.302 ; m ra 1.600 . (49) 
The ratio r s = 2Ao/8W r ° 1 is an increasing function of U/At. It changes continuously from r s — for U/At —> to 
r s = 1 for U/At — > oo. For woo < U /At < u\ it is approximately given by r s ra e -4tuo/c/ ra ther than by r s ra e ~'"^/ At / u 
for {7/4i <C 1. This is consistent with for large f /4f values (1 — r a ) being proportional to (1 — r s ) oc 4i/J7 rather than 
to (1 - r s ) oc y/At/U. 

The temperature T a of Refj^ that plays the role of our temperature T * ra 2Ao/2fc.B is plotted in Fig. 3 of that 
reference. Its [7/4i dependence is qualitatively correct. vanishes both in the limits U/At — > and [//4t — > oo. 
It goes through a maximum magnitude at an intermediate value 5/4 < U/At < 3/2. Nevertheless, the interpolation 
function used to produce it, provided in Ref. 74 of such a paper, is poor for intermediate values of U/At. However, 
that does not affect the validity of the results of Ref.—. Indeed, the studies of that reference refer to the temperature 
range < T <C T x for which the accurate dependence of T x on U/At is not needed. The goal of its Fig. 3 is merely 
illustrating qualitatively the T x behavior over the entire coupling ranged. 

The studies of RefJ^ on the square-lattice quantum liquid refer mostly to intermediate U/At values. Given the 
important role plaid by the energy scale 2Ao ra 2ks Tq in the physics of that quantum liquid, we need a more 
quantitatively accurate U/At dependence of it for intermediate U/At values. Such a dependence must be consistent with 
the qualitative physical picture of Ref.2E. According to Eqs. (gSJ) and flU) one has that 2A = 8W^ e~ As for U/At > 0. 
Here the parameter X s is a continuous decreasing function of U/At. It is given by A s = oo for U/At — > and X s = 
for U/At oo. This reveals that the magnitude X s = 1 separates two physical regimes. Consistently, A s = 1 refers to 
the U/At value U/At — w at which 2A reaches its maximum magnitude, max{2A } = [/io(uo)] e _1 = 8W® 1 (uo) e_1 - 

That as found in RefJ^ the parameter x® is for approximately uq < U/At < u\ given by x® = 2r s /n = 2e~ itu °l u /tt 
and at U/At ra u* = 1.525 reads x® ra 0.27 implies that uq ra 1.3. Moreover, in that reference it is found that 
8^ (it,) = [396.8/295] t ra 1.345 1 at U/At ra u* = 1.525. Such a result is obtained in that reference from comparison 
of the x = and m — spin excitation spectra for the high symmetry directions found by use of the c and si 
fermion description with those estimated by the controlled approximation of Ref. 3 -. On combining the equality 
Mo (u ) = SW^Uq) with the magnitude Ho(u ) ra pe 1 /^]* of Eq. 02J we find 8^ (u ) ra [2c 1 /it] t ra 1.731 1 and 
max{2A } = [n (u )] e" 1 = 8W? x (uo) e^ 1 ra 2t/?r. 

For the intermediate range {7/4i G (tio,Ui) the L//4t dependence of the energy parameter 8W^\ is of the form 
8Wgi ra [2e 1 /7r] i W(U/At). Here W(u) is an unknown function of u = [7/4t such that W(uo) = 1- On combining 
Eqs. ((48]) and ((49j one finds 2A = 8W° 1 e - 4tu °/ u . Thus fulfillment of the maximum condition d2A {u)/du = 
at u = u/At = u requires the function W[u) be given by W(u) ra (2 — u/uo) for < [(u — u )/(ui — u )} <C 1. 
In turn, the use of the above results 8^(1/*) = [396.8/295] t ra 1.345 1 and 8W°t(u ) ra pe 1 /^]* ra 1.731 1 reveals 
that the ratio 8M / s ° 1 (w*)/8H / ° 1 (u ) may be expressed as 8H / s ° 1 (m*)/8W / s 1 (mo) ra [1 — (it, — u )]. It then follows that 
W(u) ra [1 — [u — uo)] for u ra u, = 1.525. We then use a suitable interpolation function for W(u), which has 
these two limiting behaviors. The energy parameter 8W®i is for U/At intermediate values U/At £ (uq,ui) then given 
approximately by, 

8W?i ra —W(U/4t); W(u) ra 1 - (u - u Q ) e'^k ln(no) , u < [//4t < Ul , 

7T 

, 198.4 , N 

m « u* - 1 + e _1 7r— — ra 1.302; u* ra 1.525. (50 
295 

The value u ra 1.302 is that obtained from the equation 8W° 1 (u :t ) = [2e l t/ii\ [1 - (it, - u Q )\ = [396.8/295] t ra 1.731 1. 
On combining the above results we find the following approximate behaviors for the energy parameter 2Ao, 

2A ra 32te- 7T \ /I ^ U , U/At < 1 , 

= max {2A } ra 2t/Ti , U/At = u , 

ra e (1 - itUo/u) [2t/n]W(U/At), u <U/At< Ul , 

= 8W°! ra [tt [8t) 2 /2U] , U/8 2 t > 1 . (51) 
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Here W(U/At) is the interpolation function given in Eq. (|50[) . Note that d2Ao(u)/du — at u = u/At = uq, 
consistently with 2Aq reaching its maximum magnitude 2t/ir at that U/At value. The overall U/t dependence of 
Tq 2Ao/2fcs is similar to that plotted in Fig. 3 of Ref.— for T x with the U/t value at which the maximum magnitude 
is reached shifted from U/t « 5.60 to U/t s» 5.21. Moreover, that magnitude is lessened from max{T x } s» 0.625i//cs 
to max{T*} w max {2A /2/c B } w i/^fcij] « 0.318 t/fc s . 

4- Short-range incommensurate spiral spin order for 0<i<Cl and T — 

Here we provide strong evidence that for intermediate and large values of U/At and small hole concentrations 
< x <C 1 the short-range spin order of the m — ground state corresponds indeed to that of a spin-singlet 
incommensurate spiral state. In terms of the rotated-electron spins occupancy configurations the ground state is then 
a spin-singlet incommensurate spiral state for U/At > 0, m = 0, and < x -C 1. That evidence is found on combining 
several results. This includes the necessary condition for occurrence of a long-range spin order for N% — ¥ oo argued 
in Section III-C to be that the spin effective lattice is identical to the original lattice and thus N% = 2S C = N%, the 
result of Ref42 that the m = and 0<i<1 ground state is a spin-singlet state for N% 1 and thus remains so for 
N% — > oo, and the results of ReiM- concerning the spin degrees of freedom of a related quantum problem. 

As discussed in Section I, for intermediate and large values of U/At the Hubbard model on the square lattice given 
in Eqs. ([T]) and ([5]) in terms of electron and rotated-electron operators, respectively, can be mapped onto an effective 
t — J model on a square lattice with t, t' = t'(U/At), and t" = t"{U/At) transfer integrals. The role of the processes 
associated with t' — t'(U/At) and t" = t"{U/At) becomes increasingly important upon decreasing the U/At value. 
Reference 3 ^ presents rigorous results on the spin degrees of freedom of the t — J model on a square lattice with t, t', 
and t" transfer integrals. The investigations of that paper refer to small values of the hole concentration < i « 1 
and spin density m — 0. Their starting point is a suitable action first introduced in Ref 42. The use in Ref 41 of a 
staggered CP 1 representation for the spin degrees of freedom allows to resolve exactly the constraint against double 
occupancy. Within our description of the problem, this is is equivalent to performing the electron - rotated-electron 
unitary transformation. In order to achieve the rigorous result that for small hole concentrations there occurs a 
incommensurate-spiral spin order, the effective action for the spin degrees of freedom is reached after integrating out 
the charge fermionic degrees of freedom and the magnetic fast CP 1 modes. Importantly, the dependence on the hole 
concentration of the coupling constants of the effective field theory is in Ref. 34 obtained explicitly for small x. 

We consider the mapping between the above effective t — J model and the Hubbard model of Eqs. ([I} and §5§ in 
the subspace with vanishing rotated-electron double-occupancy. Accounting for such a mapping, the studies of Ref4^ 
imply that for intermediate and large values of U /At, spin-density m = 0, and small hole concentrations < i « 1 the 
ground state of the Hubbard model on the square lattice is an incommensurate spiral state. This is a rigorous result. 
However the studies of Ref4^ are not conclusive on whether for N% — > oo and 0<i<1 the m — incommensurate 
spiral ground state has short-range or long-range spin order. 

The necessary condition for occurrence in the limit N% -4 oo of a ground-state long-range spin order argued in 
Section III-C to be that the spin effective lattice is identical to the original lattice and thus N% = 2S C = N% is 
not fulfilled for < x -C 1. We thus argue that in the thermodynamic limit N% — >• oo the < i « 1 and m = 
incommensurate spiral ground state has no long-range spin order. Hence such a ground state has the same symmetry 
and basic properties both for N% ^> 1 very large but finite and for N% — >• oo. One then combines the result of Ref 42 
that for JVf 3> 1 and thus for N% —> oo the m = and < x <C 1 ground state of the Hubbard model on the 
square lattice ([T]) is a spin-singlet state with that of Ref4^ that it is a incommensurate spiral state. This consistently 
implies that for N% — > oo and intermediate and large values of U/At the m — and < x <C 1 ground state of the 
Hubbard model on the square lattice is a spin-singlet state with short-range incommensurate-spiral spin order and 
strong antiferromagnetic correlations. 

Finally, a ground-state short-range spin order does not preclude the occurrence of a ground-state long-range dimer- 
dimer order, as in a valence-bond solid. However, whether the spin degrees of freedom of the square-lattice quantum 
liquid refer for x > and m — to a valence-bond solid or a valence-bond liquid remains an open issue. 

IV. CONCLUDING REMARKS 

In this paper we considered a suitable one- and two-electron subspace in which the general operator description 
for the Hubbard model on a square lattice with ^> 1 sites introduced in Ref. 12 simplifies. When acting onto 
such a subspace the model refers to a two-component quantum liquid described in terms of charge c fermions and 
spin-neutral two-spinon si fermions. The one- and two-electron subspace can be divided into smaller subspaces that 
conserve S c and S s . Those are spanned by energy eigenstates whose generators have simple form in terms of c and 
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si fermion operators, as given in Eqs. (|28[1 and (j2T))) . When expressed in terms of c and si fermion operators, the 
Hubbard model on a square lattice in the one- and two-electron subspace is the square-lattice quantum liquid further 
studied in Rcf. 19 . That in such a subspace the c and si fermion momentum values are good quantum numbers plays a 
key role in the investigations of that reference. The one- and two-electron subspace considered in this paper contains 
nearly the whole spectral weight generated from application of one- and two-electron operators onto the exact ground 
state. 

There is a large consensus that in the thermodynamic limit N 2 — > oo long-range antiferromagnetic order sets in 
in the m — ground state of the half-filled Hubbard model on the square lattic o 11 ' 26 " — . Consistently, in this paper 
strong evidence is found that provided that in the thermodynamic limit N 2 — > oo a long-range antiferromagnetic order 
occurs in the ground state of the related isotropic spin- 1/2 Heisenberg model on the square lattice, a similar long range 
order sets in in that limit in the ground state of the half- filled Hubbard model on the square lattice for U/At > 0. Our 
results indicate that for U/At > a ground-state spontaneous symmetry breaking from 50(3) x 50(3) x U{1) for large 
but finite number of lattice sites N 2 > 1 to [U(2) x U(l)]/Z$ = [50(3) x U(l) x U(l)]/Z 2 in the thermodynamic limit 
N 2 — » oo occurs at x = and m = due to emergence of such a long-range antiferromagnetic order. Our analysis of 
the problem profits from the description of the ground states and low-lying energy eigenstates and corresponding state 
representations of the group 5*0(3) x 50(3) x U(l) in terms of occupancy configurations of three effective lattices. 
Our results indicate that the spin effective lattice being identical to the original lattice and thus N 2 = N 2 is a 
necessary condition for the occurrence of a ground-state long-range antiferromagnetic order as N 2 — > oo. The picture 
which emerges for N 2 — > oo is that of a ground state with long-range antiferromagnetic order for half filling and a 
short-range incommensurate-spiral spin order for 0<i<1 and thus N 2 < N 2 . 

An interesting issue is whether the quantum phase transition separating the x — and m = ground state from 
the < x <C 1 ground state at T — corresponds to a "deconfined" quantum critical point^Sr— . That both at x = 
and for x > the M s = 2S C spinons are confined within N s \ = S c spin-neutral two-spinon si fermions strongly 
suggests so. This would imply that the critical point is characterized by deconfined spin-1/2 spinons coupled to some 
emergent U(l) gauge field. Furthermore, in this paper the general rotated-electron description behind the c and si 
fermion operator representation was also used to provide evidence that the Mermin and Wagner Theorem may apply 
at half- filling to all values U /At > of the Hubbard model on the square lattice. 

The energy order parameters (jP — U tuaf oP of Eq. (pt3"f and 2|A| of Eq. P5|) have a different physical origin. 
The energy scale /jP can be used as order parameter of the T = 0, x = 0, and m = phase with long-range 
antiferromagnetic order. Indeed, it is proportional to the x = and m = sub-lattice magnetization vtlaf of Eq. 
(1351) . which vanishes for x > 0. In turn, the energy order parameter 2|A| is associated with the x > and m = 
short-range spin correlations. However, both /jP and 2|A| are at x = and for x > 0, respectively, identified in 
the studies of Ref. 1 - with the maximum pairing energy of the —1/2 and +1/2 spinons of a composite spin-neutral 
two-spinon si fermion. The results of that reference extend the m = and T = short-range spin order found in 
this paper for < x <C 1 to a well-defined range of hole concentrations < x < x*. According to these results, the 
x dependence 2|A| 2Ao(l — x/x^) given in Eq. (pQjf for < x <C 1 is valid for x £ (0, x*). This holds provided 
that approximately U/At £ (uq,u w ). Here > u\ where u\ w 1.6 is the U/At value at which r s = 1/2. For that 
U /At range the critical hole concentration x* equals the U/At- dependent parameter x° = 2r s /ir. The studies of Ref 
identify it with a critical hole concentration x* = x® above which there is no short-range spin order at T = 0. For 
x > x* and T — a spin disordered state without short-range order for which the energy scale 2|A| vanishes emerges. 
Consistently, 2|A| — > as < (x* — x) — > 0. The short-range incommensurate-spiral spin order discussed here for 
< x <C 1 corresponds then to a limiting case of the general short-range spin order that according to the investigations 
of RefJ£ occurs for < x < x* and approximately U/At € (uo,u n ). (That order occurs as well for U/At > u n , yet 
then the critical hole concentration x* may not be given by x° = 2r s /it.) 

As confirmed in that reference, for the square-lattice quantum liquid introduced in this paper the c and si fermions 
play the role of "quasiparticles" . There are three main differences relative to an isotropic Fermi liquid^. First, 
concerning the charge degrees of freedom, the non- interacting limit of the theory refers to At 2 /U — > rather than 
to the limit of zero interaction U 0. Second, in the At 2 /U —> limit the c fermions and si fermions become the 
holes of the "quasicharges" of Refi^ and spin-singlet two-spin configurations of the spins of such a reference rather 
than electrons. Indeed, only the charge dynamical structure factor becomes that of non-interacting spinless fermions. 
In turn, the one-electron and spin spectral distributions remain non-trivial. Third, for U/At > the si band is full 
for initial m = ground states and displays a single hole for their one-electron excited states. As found in Ref.— its 
boundary line is anisotropic, what is behind anomalous one-electron scattering properties. Those involve the inelastic 
scattering of c fermions with momenta near the isotropic c Fermi line with si fermions with momenta in the vicinity 
of the anisotropic boundary line. 

Concerning the relation to previous results on the Hubbard model on the square lattice and related models by 
other authors, as discussed above in this paper our results are consistent with and complementary to those of 
Refs i 29 i 30 ' 34 i 36 ' 38 . Elsewhere evidence is provided that upon addition of a weak three-dimensional uniaxial anisotropy 
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perturbation to the square- lattice quantum liquid, its short-range spin order coexists for N% — > oo, low temperatures, 
and a well-defined range of hole concentrations with a long-range superconducting order. 
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Appendix A: Results on the c fermion and si bond-particle description needed for our studies 

For the LWS subspace defined in Ref.— the rotated-electron operators of Eq. Q can be expressed in terms of the c 
fermion operators of real-space coordinate r,- and rotated quasi-spin operators = at + pi and <f-? — s% 3 + p% 3 and 

thus of spinon operators s-, and s% 3 and n-spinon operators pi and p% 3 . This is achieved by inverting the relations 

given in Eqs. ((6|) and (j9|) with the result, 

4t = 4e ( L , ■ -v • /V') • (l~ s ^-p^ ; ^ = ^+^)(4,c-^ , /i' J ,c). (Ai) 

The corresponding expressions of the rotated-electron annihilation operators are trivially obtained from those provided 
here. 

As given in Eq. ©, Fourier transform of the c fermion operators ft c generates the corresponding momentum c 
fermion operators ft c . Here the discrete momenta qj are good quantum numbers^. In turn, the operators ft av of 
the av fermions also introduced in Ref. 12 act onto subspaces with constant values for the set of numbers S a , N au , and 
{N av i} for v' > v and equivalently of the numbers S c , N av , and {N au i} for all v 1 ^ v. Here v = 1, C a is the number 
of 77-spinon (a = rf) and spinon {a = s) pairs confined within a composite av fermion and the maximum v value C a 
is expressed below in Eq. (|A5[) in terms of av fermion numbers N au . The above subspaces are spanned by mutually 
neutral states, that is states with constant values for the numbers of av fermions and av fermion holes. Hence such 
states can be transformed into each other by av band particle-hole processes. As a result of the transformation that 
maps the av bond particles into av fermions and provided that in the thermodynamic limit N% — > 00 the ratio 
Ncv/Nq is finite for the av branch under consideration, the latter composite fermions have long-range interactions 
associated with an effective vector potential A a ^(rj) 12 ^ 3 ^ 4 . For the one- and two-electron subspace considered in 
this paper the ratio N al ,/N 2 remains finite in the limit N 2 — > 00 only for the si fermion branch. For it the effective 
vector potential A s i(fj) reads 12 , 

7 . x \ e x „ x fry — r, ) „+ 

A s i{rj) = $o2^"*V.'i _>yT i n ^,si = ft xf^si, 

B sl {fj) = V r - x A sl (fj) = $ n ? 3 >- sl 8( fi' ~ *j) ^3 ; $o = l- (A2) 

Here B s i is the corresponding fictitious magnetic field, e X3 is the unit vector perpendicular to the plane, and we use 
units such that the fictitious magnetic flux quantum is given by $0 = 1- It follows from the form of the effective 
vector potential A s i(fj) that the present description leads to the intriguing situation where the si fermions interact 
via long-range forces while all interactions in the original Hamiltonian are onsite. 

The theory associated with the operator description introduced in RefJ^ for the Hubbard model on the square 
lattice refers to a well-defined vacuum. For hole concentrations < x < 1 and maximum spin density m = (1 — x) 
reached at a critical magnetic field H c parallel to the square-lattice plane the c fermion operators are invariant under 
the electron - rotated-electron unitary transformation. Then there is a fully polarized vacuum |0^ s ), which remains 
invariant under such a transformation. It reads, 

|0„.) = |0„; N 2 ar ) x |0 S ; N 2 J x \GS C : 2S C ) . (A3) 

Here the ry-spin SU(2) vacuum \0 V ; N 2 ^) associated with A" 2 ^ independent +1/2 77-spinons, the spin SU(2) vacuum 
|0 S ; N% s ) with N% b independent +1/2 spinons, and the c U(l) vacuum \GS C ;2S C ) with N c = 2S C c fermions remain 
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invariant under the electron - rotated-electron unitary transformation. The explicit expression of the state \GS C ', 2S, 

lg J q,c 



of Eq. (|A3I) in terms of the vacuum \GS c ;0) appearing in Eq. (|29|) is Y\~ ft\GS c ; 0). The vacuum \GS c ;0) is part 



of the electron and rotated-electron vacuum. The form of the latter vacuum is given by Eq. (|A3j) with N% = N%, 
iVf = 2S C = 0, and thus \GS C ;2S C ) replaced by \GS c ;0). Only for a m = (1 — x) fully polarized state are the 
state \GS C ;2S C ) and the corresponding N c = 2S C fermions invariant under the electron - rotated-electron unitary 
transformation for U /At > 0. 

The av momentum band of the operators ft av is associated with a well-defined av effective lattice. The number 

N% of the av band discrete momentum values exactly equals the number of sites of such an effective lattice. It is 
given by, 

Nt av = [N av + Nt\: (A4) 

where the number of unoccupied sites reads, 

JV*, = [2S„ + 2 £ iy'-v)N av .\ = [N 2 aa - Y^(v + v'- \v-v'\)N av ,]; C a =Y j vN av . (A5) 

u' — v+l v' — l u—1 

For the particular case of v = 1 the expression (IA5|) of the number of unoccupied sites of the al effective lattices 
simplifies to, 

c„ 

Ki^[N! a ~2B a }; B a = Y^N av ; a = f 1 ,s. (A6) 

t/=i 

This number equals that of al fermion holes in the al band. 

According to the results of RefJ^, all sites of the si effective lattice of x > and m = ground states are occupied 
and hence there are no unoccupied sites. In turn, the dominant contributions to the one-electron and two-electron 
excitations involve states with one and none or two unoccupied sites, respectively. For the square-lattice quantum 
liquid the expression of the related conserved number introduced in that reference simplifies to, 

Here N denotes the number of electrons. For the Hubbard model on the square lattice in the one- and two-electron 
subspace considered in this paper the number of unoccupied si effective lattice sites and thus of si fermion holes 
in the si momentum band is a good quantum number. 

As discussed in RefJ^, the 77-spinon, spinon, and c fermion description contains full information about the relative 
positions of the sites of the 77-spin and spin effective lattices in the original lattice. Hence it turns out that within the 
N% 3> 1 limit and for finite values of the hole concentration x (and electron density n — (1 — x)) the 77-spin (and spin) 
effective lattice can be represented by a square lattice with spacing a v (and a s ) given by, 



Nn- N„ 



a; a — r\ , s . (A8) 



Moreover, provided that in the thermodynamic limit N 2 — > 00 the ratio N alJ /N 2 is finite, the related av effective 
lattices can be represented by square lattices with spacing, 

ia ; N aau > 1 , (A9) 



where v = 1, C a and a = T),s. In turn, the corresponding av bands whose number of discrete momentum values 
is also given by are well defined even when is given by a finite small number, N% = 1, 2, 3, ... 

Finally, an important and useful property found in Ref. 12 is that the rjv fermions with any number v = 1, C v of 
77-spinon pairs and sv fermions with v = 2, C s spinon pairs whose energy is given by, 

e vv = 2vfi, v = 1, Cj] ; e SL , = 2v[ib H , v = 2, C s , (A10) 

remain invariant under the electron - rotated-electron unitary transformation. Here H denotes the magnitude of a 
magnetic field aligned parallel to the square-lattice plane. Such rjv fermions and sv fermions are non-interacting 
objects. Hence their energy is additive in the individual energies of the corresponding 2v 77-spinons and spinons, 
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respectively. For U /At > such quantum objects refer to the same occupancy configurations in terms of both rotated 
electrons and electrons. 
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